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A  production  system  is  modeled  as  an  activity 
network  using  a  dynamic  activity  analysis.  The 
nodes  of  the  network  represent  component  activi¬ 
ties  which  produce  final  products,  intermediate 
products  (inputs  to  other  activities) ,  or  both, 
according  to  intensity  functions  varying  with 
time.  The  arcs  of  the  network  indicate  where 
transfers  of  intermediate  product  occur.  These 
transfers  of  intermediate  product  can  be  modeled 
two  ways:  either  as  discrete  units  on  an  event 
basis,  or  as  continuous  flows.  In  the  former 
case,  a  generalization  of  critical  path  planning 
is  obtained,  whereby  variation  in  activity  inten¬ 
sities  allows  variable  resource  load  levels  and 
production  durations.  The  intensities  provide  an 
added  degree  of  flexibility  for  production  plan¬ 
ning,  smooth  loading  of  resources,  etc.  In  the 
latter  case,  linear  programming  may  be  used  to 
allocate  resources  so  as  to  maximize  throughput, 
minimize  production  costs,  and  make  determina¬ 
tions  of  capacity  and  efficiency. 
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1.  General  Dynamic  Activity  Analysis 

The  dynamic  activity  analysis  for  production  has  previously  been  set  forth 
in  references  [10]  and  [11],  in  which  the  structure  of  the  dynamic  activity 
analysis  was  presented  as  an  example  in  the  theory  of  dynamic  production  net¬ 
works.  This  chapter  is  adapted  from  reference  [10]  to  develop  the  base  for 
concepts  and  procedures  presented  in  following  sections. 

A  production  system  is  modelled  as  a  network  of  activities  which  are 
denoted  by  Alt  ...  ,  Ay.  Figure  1.1  displays  a  typical  configuration  for  an 
activity  network.  Each  node  represents  an  activity,  which  produces  final  pro¬ 
ducts,  intermediate  products  (inputs  to  other  activities),  or  both.  The  arcs  indi¬ 
cate  activity  dependencies,  ie,  where  in  the  network  transfers  of  intermediate 
product  occur.  In  the  figure,  activity  A  ]  provides  intermediate  product  input  to 
A 3,  while  A 2  serves  A3,  A 4  and  A-,,  etc.  The  double-stemmed  arrows  indicate 
activities  producing  final  products;  in  this  case,  A5,  A6  and  A8  all  produce  final 
output.  Restrictions  on  the  network  structure  appropriate  to  modeling  certain 
kinds  of  problems  will  be  discussed  in  the  following  sections. 

In  addition  to  the  input  of  intermediate  product,  activities  require  exo¬ 
genous  inputs  for  operation,  ie,  resources  not  produced  by  other  activities.  The 
operation  of  each  activity  is  allowed  to  vary  over  a  discrete  time  grid 
t-0, 1,2 . where  on  each  interval  [(— 1 , r) ,  r=l,2,3 .  inputs  (exo¬ 

genous  and  intermediate  product)  are  applied  at  uniform  rates,  and  charged  at 
time  (/ — 1).  Outputs  of  production  during  this  interval  are  forthcoming  at  time 
r,  at  which  time  they  may  be  transferred  as  intermediate  or  final  product,  or 
held  as  part  of  accumulated  inventories. 

The  operation  of  each  activity  A,  is  measured  in  terms  of  an  intensity  June- 


Figure  1.1 

Example  of  Production  Activity  Network 
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-(/),  r-0.1.2 . 

whose  value  at  time  t ,  taken  in  conjunction  with  technical  coefficients,  indicates 
production  input  during  (r,/+l)  and  activity  output  at  time  (r+l).  In  this 
fashion,  activity  intensity  is  taken  to  be  constant  on  each  interval  U.r+1).  The 
technical  coefficients  are  defined  as  follows: 

(a)  c,(r),  r=0, 1,2,  ...  ,  /'=1,  .  .  .  ,N,  where  c,(r)  is  the  amount  of  output  of 
activity  A,  at  time  t  per  unit  intensity  of  activity  A,.  If  the  output  is  a  siz¬ 
able  discrete  unit,  c,(t)  is  understood  to  be  the  fraction  of  a  unit  com¬ 
pleted  at  the  time  t  per  unit  intensity  of  operating  A:\  otherwise,  it 

-  represents  amount  (in  bulk).  In  most  applications,  these  output  coefficients 
are  constant  in  time  or  perhaps  are  direct  functions  of  cumulative  produc¬ 
tion  (to  indicate  "learning  curve"  behavior),  but  are  expressed  most  gen¬ 
erally  as  functions  of  time. 

(b)  ajk(t),  /= 0, 1,2,  ...  ,  Ar=l ,  2,  .  .  .  ,NK,  /=  1,  .  .  .  ,N,  where  a,k{t)  is  the 
amount  of  exogenous  input  k  required  at  time  t  per  unit  intensity  of 
activity  Ah  and  the  index  k  spans  all  NK  exogenous  resources  required  by 
the  production  system.  As  with  the  output  coefficients,  these  input 
coefficients  are  typically  constant  in  time  or  perhaps  direct  functions  of 
cumulative  production  (indicating  "learning  curve”  behavior),  but  are 
expressed  most  generally  as  functions  of  time. 

(c)  a,j(t),  t=0, 1,2,...,  /=1 . Ar,  y-l,  .  .  .  ,N,  where  is  the 

amount  of  intermediate  product  from  activity  required  at  time  t  per  unit 
intensity  of  operating  Ar  (In  the  case  of  discrete  intermediate  product 
transfers,  constant  transfer  coefficients  a,j  relate  the  number  of  units  of 
intermediate  product  from  activity  A.  required  per  unit  output  of  activity 
A,.  See  Chapter  2.) 
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These  coefficients  are  nonnegative  with  the  following  restrictions: 

(a)  c,(/)>0,  r=0. 1.2 . /-I . :V. 

(b)  For  each  r^O  and  each  /  =  1 ,  .  .  .  ,  V,  there  exists  a  positive  coefficient 

aik(t )  for  at  least  one  k  e  (1 . NK}. 

(c)  For  each  system  exogenous  input  A:  e  {1,  .  .  .  ,NK),  there  exists  r^O  and 

at  least  one  activity  Ar  /  e  { 1 . A’},  such  that  aik(t)  is  positive. 

The  activity  intensities  (/)  have  natural  bounds,  since  production  cannot 
be  carried  out  instantaneously.  These  bounds,  denoted  by 

r ;(/),  /- 1, _ A,  /- 0,1,2 . 


result  from  the  available  workspace  and  other  limitations  not  considered  as  exo¬ 
genous  input.  Of  course,  particular  allocations  of  exogenous  inputs  may  limit 
intensity  further. 

A  production  plan  is  a  specification  over  some  finite  period  [0,T)  of  the 
activity  intensities 


f-0 . T- 1 

i-l . N 


Such  a  specification  determines  the  production  system  demand  histories  for 
exogenous  inputs  (resources)  as  well  as  the  output  histories.  For  example,  the 
demand  for  exogenous  input  k  during  the  interval  (r,r+l)  is  given  by 

(')*/*( ') ; 

/"I 


and  for  the  ca^e  that  activity  A  v  is  the  only  activity  producing  final  output,  and 
Ay  produces  only  final  product,  the  production  system  cumulative  output  dur¬ 
ing  [0.T)  is  given  by 


^  Z\(t)c\(t  + 1)  . 

/- 0 

To  be  feasible,  a  production  plan  must  of  course  satisfy  constraints  con¬ 
cerning  intermediate  product  transfers,  available  exogenous  inputs,  etc.,  which 
are  discussed  in  following  sections.  Expression  of  such  constraints  and  further 
development  of  the  model  depends  considerably  on  whether  the  transfers  of 
intermediate  product  are  considered  to  be  large,  discrete  units  transferred  on  an 
event  basis,  or  are  viewed  as  continuous  flows.  The  former  point  of  view  seems 
appropriate  for  detailed  modeling  of  large-scale  construction  such  as  shipbuild¬ 
ing,  in  which  virtually  no  work  on  the  installation  of  various  systems  can  begin 
until  large  assemblies  such  as  the  hull  structure  are  completed.  On  the  other 
hand,  the  latter  point  of  view  is  more  appropriate  for  production  processes 
where  either  output  has  no  discrete  nature  or  else  a  considerable  number  of 
discrete  output  units  may  be  produced  in  a  single  time  interval.  Such  would  be 
the  case  for  aggregate  plants,  food  processing  plants,  paper  production,  and 
many  other  examples. 

In  the  sections  that  follow,  the  model  is  developed  first  for  discrete 
transfers  of  intermediate  product.  Starting  with  simple  extensions  of  Critical 
Path  Methods,  the  model  is  generalized  into  a  fully  dynamic  structure  with 
intensity  histories,  output  streams,  and  product  inventories.  A  transition  is  then 
made  to  continuous  product  flows  for  the  final  development  of  the  model. 


2.  Dynamic  Production  Planning  with  Discrete  Transfers  of  Intermediate 
Products 

2.1.  Nature  of  Discrete  Product  Transfers 

For  the  discrete  model  of  product  transfers,  we  shall  assume  output  of 
each  activity  is  measured  in  discrete  units,  and  that  product  transfers  can  only 
occur  in  integral  amounts  of  such  units.  For  each  activity  An  product  transfers 
from  every  activity  Aj  supplying  A,  are  required  for  A,  to  initiate  production  of 
each  output  unit,  according  to  product  transfer  coefficients 

a,j,  i=l,  ...  ,N,  jj*i, 

w'here  a tj  is  the  number  of  output  units  from  activity  Aj  required  in  the  produc¬ 
tion  of  one  output  unit  by  Ar 

Note  that  in  the  discrete  model,  product  transfer  coefficients  are  not 
transfers  per  unit  intensity.  However,  as  the  intensity  of  activity  A ,  is  increased, 
more  frequent  product  transfers  would  be  required  of  Ar  implying  greater 
intensity  for  Aj.  Such  activity  interaction  is  entirely  similar  in  the  case  that  pro¬ 
duct  transfers  are  modeled  directly  proportional  to  activity  intensity  (continuous 
flow  model,  Chapter  3). 

Discrete  product  transfers  are  intimately  related  to  the  notion  of  pre¬ 
cedence  in  critical  path  methods.  That  is,  the  work  package  of  A,  producing  one 
unit  is  preceded  by  the  work  package  of  A}  producing  atj  units.  In  order  to 
more  fully  develop  this  relationship,  it  will  be  instructive  to  turn  to  critical  path 
techniques,  and  to  consider  extensions  of  such  techniques  in  the  spirit  of  the 
dynamic  activity  analysis. 
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2.2.  Activity  Analysis  Extensions  of  Critical  Path  Techniques 

2.2.1.  Review  of  CP\1 

PERT  and  CPM  are  network  planning  techniques  for  project  management. 
These  techniques  are  usually  applied  to  one-time  efforts,  in  which  similar  work 
may  have  been  done  previously,  but  it  is  not  being  repeated  on  a  production 
basis. 

Critical  path  methods  involve  both  a  graphical  portrayal  of  the  interrela¬ 
tionships  among  the  activities  of  a  project,  and  an  arithmetic  procedure  which 
identifies  the  contribution  of  each  activity  to  the  overall  schedule,  as  briefly  dis¬ 
cussed  below'. 

In  a  CPM  project  network,  nodes  denote  activities  which  require  time, 
manpower,  and  facilities  to  complete,  as  illustrated  in  Figure  2.1.*  The  arrows 
indicate  work  flow  dependencies,  eg,  activity  A2  must  be  completed  before 
activity  A4  is  performed.  To  maintain  correct  network  logic,  cycles  are  prohi¬ 
bited.  An  activity  is  thus  any  portion  of  the  work  that  may  not  begin  until  other 
portions  are  completed. 

Each  CPM  activity  A,  has  a  time  assignment  t,  which  is  an  estimate  of  the 
duration  required  to  complete  that  portion  of  project  work  represented  by  the 
activity.  The  basic  scheduling  computations  are  a  forward  pass  and  a  backward 
pass  through  the  network.  Typically,  start  times  for  the  initial  activities  (in  the 
figure,  Ax  and  A2)  are  given  data,  so  that  the  forward  pass  is  performed  first. 
Based  on  the  given  start  times  and  activity  durations,  the  forward  pass  com¬ 
putes  the  earliest  start  and  earliest  finish  times  for  each  node  (activity).  Similarly, 
when  finishing  times  for  final  nodes  (in  the  figure.  A7  and  .*f3)  are  specified  or 


*  In  most  presentations  of  CPM.  the  roies  of  nodes  and  arcs  are  recersed  from  the  roles  -seo  here  Either 
formal  can  he  used  in  practice.  See  reference  [5!. 


computed  by  the  forward  pass,  the  backward  pass  can  compute  the  latest  start 
and  latest  finish  times  for  each  node.  The  times  are  latest  in  the  sense  that  if 
any  were  later,  project  finishing  times  would  not  be  met. 

Comparison  of  backward  and  forward  pass  results  will  indicate  that  certain 
paths  through  the  network  will  be  critical.  That  is,  start  and  finish  limes  for 
activities  on  these  paths  any  later  than  the  earliest  times  will  cause  later  project 
finishing  times.  Other  paths  will  have  slack,  whereby  nodes  on  such  paths 
could  have  larger  time  assignments  without  impacting  these  finishing  times. 
Alternatively,  the  starting  times  for  such  nodes  (activities)  could  be  delayed. 

In  general,  the  critical  path  P  between  an  initial  node  and  a  final  node  is 
the  longest  path  through  the  network  connecting  those  nodes,  where  the  length 
of  the  path  T(P)  is  the  sum  of  the  time  assignments  of  nodes  on  the  path  P. 
The  slack  of  an  alternate  path  P'  between  the  same  two  nodes  is  given  by 
T(P)—T(P  ).  The  scheduling  slack  s,  of  a  node  A,  is  the  minimum  slack  value 
of  all  slack  paths  containing  that  node,  and  is  given  by  the  difference  between 
the  late  and  early  start  (or  finish)  times  for  Ar 

Associated  with  each  activity  are  requirements  for  manpower  and  other 
resources.  A  determination  of  start  and  finish  times  for  all  the  activities  allows 
the  computation  of  resource  demands  by  the  project  for  a  given  schedule.  The 
variation  in  resource  loading  resulting  from  shifting  slack  activities  between 
early  and  late  start  can  be  studied.  See  references  [5,6.12], 

2.2.2.  Activity  Intensities 

The  fixed  time  and  resource  assignments  to  each  activity  in  critical  path 
techniques  imply  fixed  application  rates  of  resources.  Ordinarily,  in  a  construc¬ 
tion  project,  one  has  the  flexibility  to  vary  such  rates  to  accomplish  leveling  of 
resource  utilization,  balancing  of  interdependent  processes,  etc.  As  a  first 
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extension  of  critical  path  methods,  we  shall  make  the  resource  assignments  to 
each  activity  variable,  according  to  the  following  restricted  dynamic  activity 
analysis. 

We  consider  a  critical  path  network  of  N  activities  . Av  modeling  a 

construction  project  in  which  A rK  exogenous  resources  are  utilized.  The  rate  of 
application  of  resources  to  each  activity  A,  can  be  modeled  linearly  in  terms  of 
a  strictly  positive  intensity  variable  r,  taken  with  technical  coefficients  defined  as 
follows: 

(a)  c,  =  fraction  of  activity  A,  completed  per  unit  intensity,  per  unit  time, 
/'- 1,  .  .  .  ,N. 

(b)  alk  =  amount  of  exogenous  resource  k  applied  to  A ,  per  unit  intensity,  per 

unit  time,  /=1 . N,  k= 1 . NK. 

The  intensity  assignment  z,  to  activity  A ,  is  maintained  at  a  constant  value 
between  start  and  finish  times  for  A,,  but  the  particular  level  of  intensity  is  a 
decision  variable.  In  this  way,  the  time  assignment  to  each  activity  A,  becomes 
a  variable,  given  by 


Typically,  for  each  activity  A,,  there  is  some  maximum  practical  rate  at 
which  that  activity  may  operate,  implying  an  intensity  bound  :r  These  bounds 
reflect  the  conditions  peculiar  to  each  activity  such  as  available  workspace. 

Up  to  this  point,  this  analysis  resembles  CPM  time-cost  trade-off  pro¬ 
cedures  (see  reference  [5]),  in  which  activity  durations  may  vary  between 
crash  and  normal  durations,  and  there  is  a  linear  relation  between  activity 
duration  and  costs  of  resources  used.  However,  in  our  formulation,  the  rela¬ 
tionship  between  the  level  of  resources  demanded  by  each  activity  and  activity 


duration  is  made  explicit,  so  that  between  the  start  and  finish  times  for  an 
activity  the  rate  of  application  of  resource  k  to  activity  A,  is  given  by 


z,atk  =  — ,  A  =  1 . NK,  /= 1 . /V.  (2.2) 

tic, 

Pictorially,  the  "resource  box"  assigned  to  an  activity  can  be  "flattened  and 
lengthened"  (while  keeping  the  area  constant)  by  decreasing  activity  intensity, 
as  shown  in  Figure  2.2. 

By  developing  technical  coefficient  data  and  defining  intensities,  one  can 
then  vary  activity  intensities  to  accomplish  time  substitution  of  activity  loading. 
This  approach  contrasts  with  pure  scheduling  strategies,  which  seek  to  minimize 
resource  peaks  by  shifting  a  fixed  box  for  each  activity  between  early  and  late 
starting  times.  (See  references  [3]  and  [5]  for  discussions  of  scheduling  stra¬ 
tegies;  see  sections  2.6  and  2.7  for  investigations  of  smoothing  and  resource 
leveling.) 

An  initial  selection  of  activity  intensities 


allows  the  ordinary  CPM  scheduling  calculations  of  early  (late)  start  and  finish 
times  and  activity  slack,  using  the  node  durations 


In  particular,  operating  all  activities  at  their  intensity  bounds,  ie. 

-0  -  />i  v 

•/  1  1  * . . 


yields  a  minimum  project  time  span. 
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The  intensity  -°  of  activity  A,  is  said  to  be  critical  if  the  lime  assignment  to 

A,  , 


is  such  that  activity  A,  is  on  a  critical  path  in  the  activity  network.  That  is,  a 
lower  intensity  assignment  to  activity  A ,  would  extend  the  project  time  span. 

Any  activity  A ,  may  be  made  critical  by  lowering  its  intensity  assignment 
o  to 


where  s,°  is  the  slack  of  activity  A ,  calculated  by  the  CPM  scheduling  computa¬ 
tions  for  the  time  assignments  (2.4),  Other  activities  on  the  same  slack  path 
would  then  also  be  made  critical,  as  their  slack  would  be  eliminated.  In  general, 
the  slack  of  a  given  slack  path  through  the  network  can  be  allocated  among  the 
activities  on  the  path,  whereby  the  reduced  intensity  assignments  would 
become  critical.  Eliminating  all  activity  slack  in  the  network  leaves  all  activities 
critical. 

Reducing  activity  intensities  to  criticality  tends  to  smooth  project  resource 
loading  histories,  as  the  resource  demands  of  slack  activities  are  reduced  to 
lower  levels.  See  the  discussion  and  examples  in  sections  2.6,  Smooth  Loading, 
and  2.7,  Resource  Leveling.  See  also  reference  [11]  for  further  examples. 
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2.3.  Extension  to  Output  Streams 

Critical  path  techniques  are  oriented  to  the  single  construction  project 
scenario.  But  consider  a  production  system  producing  a  stream  of  outputs,  such 
as  a  shipyard  producing  a  sequence  of  ships.  One  could  develop  an  activity  net¬ 
work  for  each  ship,  plan  the  resource  loads,  and  then  attempt  to  overlay  these 
loads  to  predict  overall  yard  resource  demands. 

But  desirable  elements  of  the  work  flow  and  planning  flexibility  may  be 
missing  from  this  approach.  Consider  the  simplistic  network  shown  in  Figure 
2.3  for  producing  a  ship  (activities  are  nodes,  arcs  show  precedence).  After 
activity  A  finishes  work  on  the  first  ship,  activity  B  and  activity  C  start  work  on 
the  first  ship.  But  in  view  of  the  sequence  of  ships  to  be  produced,  activity  A 
w'ould  then  commence  work  on  the  second  ship.  An  extended  network  is 
immediately  suggested,  whose  nodes  may  be  denoted  by  (/,/»),  w’here  /  denotes 
the  activity,  and  m  denotes  the  output  unit.  For  three  ships  produced  by  the 
above  network,  one  has  the  extended  netw-ork  shown  in  Figure  2.4. 

There  can  now  be  an  intensity  assignment  ztm  to  activity  A,  working  on 
output  unit  m,  which  corresponds  to  a  time  assignment 


hm 


to  node  (i,m)  of  the  extended  network. 

Technical  coefficients  can  also  vary  by  output  unit  to  admit  learning  curve 
behavior,  whereby  either  the  level  of  resources  required  per  unit  intensity 
declines  or  the  output  per  unit  intensity  per  unit  time  increases.  The  coefficients 
are  modified  as  follows: 


IL 
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(a)  cim  ~  faction  of  output  unit  /?;,  activity  A,  completed  per  unit  intensity, 
per  unit  time: 

(b)  a>km  =  amount  of  resource  k  applied  to  activity  A,  working  on  output  unit 
m  per  unit  intensity,  per  unit  time. 

The  time  assignment  to  activity  A,,  output  unit  /??,  is  then  given  by 

'im  -  >  (2.1) 

- 1  m  ctm 

and  the  rate  of  application  of  resource  k  is  given  by 

-im^iknr 

The  extended  network  concept  allows  the  integrated  planning  of  a  series  of 
construction  projects,  effectively  treating  such  a  series  as  a  single  project. 
Rather  than  attempting  to  overlay  or  to  sequence  independently-treated  CPM 
networks,  one  could  use  an  extended  network  for  overall  resource  management 
and  allocation.  See  section  2.7  for  an  example  of  resource  leveling  on  a  two- 
ship  extended  network. 
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2.4.  General  Dynamic  Model  with  Discrete  Transfers 

We  now  return  to  the  dynamic  activity  analysis,  whereby  activities  have 
intensity  time  histories 


defined  over  a  discrete  time  grid  r=0, 1,2,  •  •  •  .  Technical  coefficients  and 
intensity  bounds  are  also  allowed  to  be  variable  with  time.  We  shall  assume 
that  activity  output  units  are  so  large  that 

z,U)c,(t)  <  1  (2.8) 

for  all  activities  A,,  i—  1,  .  .  .  ,  N,  and  t— 0, 1,2,  •  •  •  .  (If  such  is  not  the  case, 
a  finer  time  grid  may  be  considered.) 

We  shall  also  assume  the  activity  network  representation  of  the  production 

system  is  cycle-free.  Product  transfer  coefficients  aIJ  are  taken  to  be  constant, 

so  that  the  arcs  of  the  network,  while  identifying  intermediate  product 
transfers,  also  indicate  the  precedence  of  work  flow.  An  arc  from  activity  A;  to 
activity  A ,  means  that  that  Ay  must  complete  (n)dlJ  output  units  before  activity 
A,  can  start  production  of  its  mh  output  unit. 

By  suitable  modification,  the  extended  network  discussed  in  section  2.3  can 
be  used  to  explicitly  display  each  discrete  product  transfer.  For  each  activity  .4 
and  output  unit  m,  the  immediate  predecessors  of  node  (i,m)  in  the  extended 
network  are  revised  to  be 
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—  1),  except  when  m  =  1;  and 

(J.a.jHi),  except  when  <7((=0. 

In  section  2.3,  the  extension  of  a  critical  path  network  in  order  to  model  output 
streams  is  an  example  of  an  extended  network  in  which  all  a,(  equal  zero  or 
one. 

Arcs  in  the  activity  network  are  allowed  to  represent  both  critical  path-style 
precedences  or  true  product  flow.  For  example,  suppose  activities  A ,  and  At 
occupy  the  same  workspace,  but  the  installation  of  A ,  product  is  required 
before  the  installation  of  Aj  product  may  occur.  In  terms  of  the  model,  we  can 
think  of  activity  Aj  requiring  the  workspace  with  activity  A ,  product  installed  as 
intermediate  product  input.  With  this  interpretation,  an  activity,  which  has  a 
critical  path-style  precedence  relationship  to  m  successor  activities,  simultane¬ 
ously  produces  m  distinct  products  according  to  its  intensity.  Each  of  the  m  arcs 
emanating  from  the  activity  represents  the  transfer  of  a  distinct  product. 

For  other  activities,  the  conception  of  product  flow  may  be  more  immedi¬ 
ate,  such  as  the  case  where  activity  A ,  produces  subassemblies  to  be  integrated 
into  the  structures  produced  by  activity  Aj.  It  may  be  the  case  that  the 
subassemblies  produced  by  A,  supply  several  activities,  so  that  arcs  emanating 
from  A,  show  flows  of  the  same  product  to  be  shared  by  successor  activities. 
Calculations  for  each  of  these  cases  of  product  How  are  presented  in  section  2.5. 

Inventories  of  intermediate  products  may  be  subject  to  storage  capacities. 
In  general,  there  is  a  storage  capacity  for  each  product  type.  For  the  presenta¬ 
tion  of  calculations  in  the  next  section,  we  shall  initially  assume  that  the  arcs  of 
the  network  represent  critical  path-style  precedences,  ie,  each  arc  represents  the 
transfer  of  a  distinct  product.  Accordingly,  we  denote 

capii.j),  /  =  ! . .V,  j  >  /, 
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to  be  the  capacity  of  stored  output  of  activity  At  to  be  used  by  activity  A.. 
where  stored  output  is  the  number  of  units  started  by  A,,  less  the  number  of 
output  units  started  by  A,  times  ajr*  To  be  feasible,  a  production  plan 


| /— 0 . T- 1 


must  satisfy  inventory  capacity  constraints  as  well  as  constraints  governing  the 
adequacy  of  intermediate  product  transfers.  These  are  discussed  section  2.5. 


■  In  the  case  that  capacity  of  product  storage  is  unconstrained,  the  capacity  may  be  set  equal  to  some  '-ery 
large  number  for  programming  purposes. 
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2.5.  Dynamic  Production  Correspondences 

Dynamic  production  correspondences  relate  time  histories  of  final  outputs 
to  time  histories  of  exogenous  resource  demands.  Their  abstract  structure  is 
discussed  in  references  [7,8,9], 

Two  kindsof  correspondence  calculations  appropriate  for  production  plan¬ 
ning  are  investigated  in  this  section.  The  presentation  is  a  generalization  of 
that  in  reference  [11],  in  which  product  transfer  coefficients  were  taken  to  be 
binary,  and  product  inventories  were  not  explicitly  considered. 

2.5.1.  Forward  Correspondence 

The  forward  correspondence  addresses  the  following  question:  for  given 
time  histories  of  the  input  of  exogenous  resources,  determine  the  time  histories 
of  output  unit  completions  possible.  Let 

?(/)-■*!(*) . XyK(i)  ,  r=0 . 7-1,  (2.9) 

be  the  given  time  histories  of  resources  available  for  production  input  during 
[0,7).  We  consider  a  preallocation  of  each  resource  k  among  the  activities, 
denoted  by 

XO(t)  =  XO^t) . Xdxf;(t)  ,  t=0 . 7-1,  (2.10) 

where 

XOk(t)  =  XOlk(t ) . XOyk(t)  .  k= 1 . A IK, 


t= 0 . 7-1, 


(2.11) 
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and 

Y.XO.Jt)  =  .rA (/)  .  A=1 . XK,  t= 0 . T-l.  (2.12) 

,  =  i 

The  full  possibilities  of  production  result  from  the  union  of  the  sets  of  network 
output  histories  obtainable  from  all  preallocations  XO(t)  of  T(t). 

An  A-by-.Y  incidence  matrix 

D(i.j) 

will  be  used  to  describe  the  network,  with  D{i,j)  =  \  if  there  is  an  arc  in  the 
activity  network  from  activity  A,  to  activity  Af  (ie,  dJt  >  0),  and  D(/.j)  = 0  if 
not.  As  the  network  is  cycle-free,  we  shall  assume  the  activity  ordering  is  such 
that  A,  may  serve  only  activities  Aj.  j>i,  /  =  1 . X. 

The  following  notation  will  be  convenient.  Let 
Afr(i)  =  number  of  output  units  completed  by  activity  A ,  by  time  / ; 
and  X,(i)  =  number  of  output  units  started  by  activity  A.  by  time  t. 

With  this  notation,  when  activities  A,  and  Aj  are  connected  by  an  arc  (ie. 
D(i.j)  =  1),  the  time  history  of  intermediate  product  inventory  of  A,  to  be  used 
by  Aj  is  given  by 

in v, ( i.j )  =  M, ( / )  -  a j, X, <j).  (2.13) 

The  constraints  on  intensity  functions  for  the  forward  correspondence  are 
as  follows: 

FdL  :,(t)  ^  0.  /=1 . X.  /=  0 . 7—1. 


F(2).  :,(r)  <  :,(t).  /-l 


V,  t=  0 


T-l. 
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F<3).  fl/A. (/)-(/)  ^  A'0,A(r).  /- 1 . V,  /=0 . T- 1. 

F(4).  j„c\(/  +  l)-(/)  ^  M,(j)  -  I^c/7+l)-(r), 

r=*0 

/  =  1 . N,  all  j  such  that  D(J,i)=l,  /=1 . T—\,  and 

a,jC,(  l)r,(0)  <  0,  '/  =  1,  .  .  .  ,  X,  all  j  such  that  D(j,i)  =  \. 

F(5).  c,(/+ \)z,(t)  <  capii.l)  +  ahN,+ j(/)  -  £  c,(r+l)  r,  (r), 

7=0 

/  =  1 . /V,  a//  /  such  that  D(i,l)—\,  t= 1,  ....  7"—  1,  and 

c,(l)z,(0)  ^  cap(i,l),  /=1,  .  .  .  ,  X,  all  I  such  that  D{i,D  —  1, 
where 


M,{i)  =  .Max 


n  integer |  /? <  £  c((r)j((--l) 

7=1 


/=i . r. 


Mq(i')  =  0, 


and 


;V,  ( / )  =  M/7 


/7  integer 1  77  ^  C  1 ) 

7=1 


f=i . r, 


(2.14) 

(2.15) 


(2.16) 


/V0(/)  =  0. 


(2.17) 
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The  constraints  F(l)  and  F(2)  allow  non-negative  values  for  intensities 
less  than  or  equal  to  the  upper  bounds.  Constraints  F(3)  limit  intensities  to 
levels  allowed  by  preallocation.  Constraints  F (4)  limit  intensities  to  levels  al¬ 
lowed  by  available  intermediate  products,  while  constraints  F(5)  insure  that  in¬ 
termediate  product  inventories  do  not  exceed  capacity .- 

A  production  plan  satisfying  constraints  F(l)  through  F<5)  is  said  to  be 
forward  feasible.  Clearly,  the  output  histories  one  can  obtain  are  not  unique. 
As  an  initial  approach,  we  shall  consider  a  greedy  policy,  ie,  when  an  activity  is 
called  upon  to  operate,  it  runs  at  maximal  intensity.  Such  a  policy  may  be  cal¬ 
culated  forward  in  time,  performing  the  four  steps  of  the  algorithm  shown 
below  at  each  instant  of  time.  It  will  be  convenient  to  keep  track  of  cumulative 
production  using  the  variables 

camfi),  /  =  !,  .  .  .  ,N,  /  =  1,2,  .  .  .  ,T. 


where 


cumt(i)  =  £:,(;)r,(7+I).  (2.18) 

r-0 

The  greedy  algorithm  for  the  forward  correspondence  is  then  as  follows: 

Step  FO.  Initialize  M^i)  =  ,Vq(/')  =  citm0(i )  =  0.  /  =  1 . V, 

and  set  t=0. 

Step  FI.  For y'=l . N,  setZj(t)  = 


\ 


Min. 


:,(t),  Min 


A-!, 


YA 


XOJk(t) 


a,, it) 


Min 

all  i 

such  that 
D  <  /,./ )  “  1 


Ml(i)-allcum,(j) 
aJlc/(t+ 1 ) 


Step  F2.  Set  Z\(t)  =  iy(/). 


Step  F3.  /"or  i=iV—\,N—2,  ....  1.  set  :,(()  = 


Min 


zit ),  Min 

all  i 

such  that 
D  ( i.j )  - 1 


capii.j )  +  <?y,|/?(/)+Ar,(/>]  -  cum,(i ) 


c,(/+l) 


where 


£(/) 


1  if  r/(/)c/-(r+l)+CM/M,0‘)  >  Ar,(/) 
0  if 


(2.19) 


an  indiciator  variable,  shows  whether  activity  Aj  initiated  production  of  a  new 
output  unit  during  the  current  time  period. 


Step  F4.  Set  =  N,(j)  +  /3 O’),  oil  j, 


where  P(j)  is  defined  by  (2.19).  Set 

M,+\(J)  =  M,(j)  +  a  O’),  oil  j, 


where 


a(j) 


1  if  Zj(t)Cj(t+ 1)  +  cum,(j)  >  A/,(/)+ 1 
0  if  not 


(2.20) 


is  an  indicator  variable,  showing  whether  activity  completed  production  of 
an  output  unit  during  the  current  time  period.  Set 


cuml+l(i)  =  cumt(i)+z,(t)Cj(i+ 1),  /  =  1 .  .  .  .  ,  .V. 


Increment  r-— r+1.  Stop  if  i=  T\  otherwise,  go  to  Step  1. 

Steps  FI  and  F2  insure  the  intensities  do  not  exceed  levels  allowed  by 
intensity  bounds,  available  exogenous  inputs  and  available  intermediate  pro¬ 
ducts.  Step  F3  insures  that  intermediate  product  inventories  do  not  exceed 
capacity.  Step  F4  updates  inventories  and  increments  time. 

From  the  calculated  intensity  histories,  one  can  observe  the  earliest  start¬ 
ing  and  finishing  times  of  each  output  unit  by  each  activity.  In  fact,  the  array 


cum,  {t) 


(-1 . v 

r-0 . T 


indicates  the  time  unit  by  time  unit  progress  of  each  activity. 

One  can  also  obtain  the  system  use  of  resources  from  the  intensity  his¬ 
tories.  During  the  time  period  lr,r-fl),  the  actual  use  of  resource  A' is  given  by 

l=,(')alkU). 

/-I 


In  the  case  where  activity  output  supplies  both  intermediate  and  final  (net) 
product,  or  in  the  case  where  activity  product  is  shared  by  two  or  more  follow- 
on  activities,  a  preallocation  must  also  be  made  of  such  output.  The  constraints 
on  intensity  functions  and  the  forward  greedy  algorithm  can  easily  be  modified 
to  handle  the  above  cases. 

As  an  example,  suppose  activity  A-  produces  one  product  which  supplies 

final  output  as  well  as  activities  Ay  and  Ay+l.  Let  5 _(/),  /— 1, 2 . T, 

denote  the  fraction  of  the  output  history  of  A _  going  to  final  output  at  time  t. 

where  0^5.,.(r)^l.  Let  A_((r),  /‘=T(2 . T,  i**y,  y+1,  be  non-negative 

real  numbers  denoting  a  preallocation  of  A.  output  at  time  t  among  v  and 
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as  intermediate  product  input.  These  coefficients  satisfy 

±vy(i)  +  Affy+1(/)  =  1  -  5_(/),  r-1 . T.  (2.21) 

For  simplicity  of  exposition,  we  shall  assume  no  activities  other  than  A _  supply 
Ay  and  A  +l.  We  denote 

cap(  tt) 


to  be  the  capacity  of  stored  output  of  .4_  to  be  used  by  Ay  and  Ay^x,  where 
stored  output  is  the  number  of  output  units  for  intermediate  uses  started  by  An 
less  the  number  of  output  units  started  by  Ay  times  ay7T,  and  less  the  number 
of  output  units  started  by  Ay+X  times  ayJrX  n. 

The  forward  correspondence  constraints  F(l),  F(2),  and  F (3)  on  intensity 
functions  require  no  modification,  but  constraints  F(4)  and  F(5)  must  be 
altered  as  follows: 


F(4)'.  Same  as  F(4)  for  i^y,  y+1.  For  i=y,y+ 1,  we  have 

a^C'it+DZ'it)  <  A/, (it,/)  -  X  a(iTf((r+l)'(r), 

r-0 

r=l . T—  1,  and 

a/77c,(l)r,(0)  <  0. 

F (5) Same  as  F(5)  for  i^rr.  For  i=i t,  we  have 

l!  ^/(f+DcJr+DzJ/)  ^ 

l-y 

cap(TT)  +  X  (/)  -  I  I  A./lr+Dc-fr+l):,,!*), 

/" y  l~y  r— 0 
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Step  F2  .  Same  as  step  F2. 


Step  F3'.  Same  as  step  F3.,  except  for  /=-,  we  have  zv(t)  = 

Minimum  of  f_ ( / )  and 

cap{rr)  +  ^  ajn  [/30)+.V,0)l  ~  I'll^/r+Dc^r+DrJr) 
_ j-y  1 _ J  ,/- yr-0  _ 

X  A^/r+Dc^Cr-f  1) 
j“y 

where  (3(J)  is  defined  by  (2.19). 

Step  F4’.  Same  as  step  F4.,  except  in  place  of  we 

set 


.  Ml+](7r,l)  =  +  a(~,l),  l=y,  y+\, 


where 


oc{tt,I)  = 


1  if  X  An./(r+l)cff(r+l)rjr(7)  >  A/,(7r,/)+l 

T-0 

0  if  not. 


(2.26) 


Production  systems  in  which  an  activity  produces  several  products,  each  of 
which  serves  several  follow-on  activities,  are  handled  in  an  analagous  manner, 
whereby  each  product  must  be  preallocated  among  uses.  For  simplicity  of 
exposition,  in  algebraic  presentations  in  the  remainder  of  Chapter  2,  it  is 
assumed  that  each  arc  in  the  production  network  represents  the  transfer  of  a 
distinct  product  (eg,  critical  path-style  precedences). 

The  greedy  policy  can  also  be  calculated  starting  at  a  time  point  with  pro¬ 
duction  in  progress.  This  is  done  by  appropriately  initializing  the  arrays  of 
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cumulative  activity  production. 


cum,(i) 


M,(i) 


and 


in  step  FO.  of  the  algorithm.  In  this  way  production  replanning  required 
because  of  work  interruptions  or  design  changes  can  be  accomplished. 

As  the  forward  greedy  policy  assigns  activity  intensities  at  the  upper 
bounds  allowed  by  constraints  F(l)  through  F(5)  at  each  instant  of  time,  the 
throughput  of  production  is  maximal  for  the  preallocations  made.  The  full 
spectrum  of  output  histories  obtainable  from  the  vector  X(t)  of  histories  of 
exogenous  inputs  is  only  realized  when,  in  addition  to  the  alternatives  for  preal¬ 
locating  exogenous  input  histories  and  transfers  of  activity  outputs  as  inter¬ 
mediate  products  and  net  outputs,  one  considers  also  the  operation  of  activities 
with  intensities  less  than-the  maximal  values  calculated  by  the  greedy  policy.  In 
particular,  the  preallocations,  in  conjunction  with  the  intensity  bounds,  will 
ordinarily  be  imperfectly  balanced  so  that  some  activity  intensities  may  be 
operated  at  less  than  the  maximal  intensities  during  certain  periods  and  still 
yield  the  net  output  streams  calculated  by  the  greedy  policy.  Such  shifts  in  the 
intensity  functions  and  the  related  input  histories  represent  time  substitutions  of 
the  input  histories  applied.  Time  substitutions  are  studied  in  sections  2.6, 
Smooth  Loading,  and  2.7,  Resource  Leveling. 

In  the  special  case  that  product  inventories  are  not  capacity-constrained, 
and  intensity  bounds,  technical  coefficients,  and  resource  preallocations  are  con¬ 
stant  in  time,  the  forward  greedy  policy  is  equivalent  to  a  CPM  forward  pass 
through  the  extended  network,  with  time  assignments 


c 


Mm 


Min 
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xo , 


a,k 


(2.27) 
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to  each  node  (/>»)  of  the  extended  network.  This  result  also  holds  when 
technical  coefficients  and  intensity  bounds  vary  by  output  unit,  as  discussed  in 
section  2.3. 

2.5.2.  Backward  Correspondence 

The  backward  correspondence  addresses  the  following  question:  determine 
the  activity  intensities  so  that  given  dated  sequences  of  net  output  units  may  be 
obtained.  We  shall  assume  one  or  more  activities  produce  final  outputs. 
Retaining  the  property  that  any  activity  A,  can  supply  only  activities  Aj  for  j  >  i, 
the  ordering  of  the  activities  may  be  made  so  that 
A,,  /'— N,N—  1,  .  .  .  ,  N-p+1  (p^l)  yield  final  (net)  outputs. 

In  the  case  two  or  more  such  activities  produce  the  same  final  output,  we 
shall  assume  a  preallocation  of  net.  output  requirements  is  made  among  these 
activities,  so  that  one  may  express  output  requirements  as  time  histories 

U((t) 

of  cumulative  output  that  must  be  completed  by  each  final  activity 
A ,,  i=N-p+ 1.  .  .  .  ,Ar.  For  example,  suppose  activity  Ay  must  complete  two 
units  by  time  25,  one  more  unit  by  time  30,  and  a  fourth  unit  by  time  45. 

Then  the  time  history  j^Tr)  •  is  constructed  as  follows: 

4  if/^45 
3  if 30<r<45 
2  if  25 <30 
0  if  r<25  . 


(2.28) 
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The  significance  of 


l\(r) 


is  clear;  if  the 


at  a  given  time  t  does  not  equal  or  exceed 
infeasible. 


cumulative  production  of  activity 
U\(t ),  the  production  plan  would  be 


As  a  first  approach  it  will  also  be  assumed  that  none  of  the  activities  pro¬ 
duce  the  same  intermediate  product,  ie,  there  are  no  alternative  processes  in 
the  production  network.  This  simplifies  the  calculations,  because  then  one 
need  not  attempt  to  apportion  among  alternative  activities  production  of  pro¬ 
duct  needed  by  follow-on  activities. 

Calculations  for  the  backward  correspondence  are  made  backwards  in  time 
from  some  time  horizon  T,  by  which  time  all  output  must  be  completed.  We 
again  assume  a  preallocation  of  the  services  of  fixed  capital  equipment  and 
other  resources,  denoted  by 


XO,(t) 


k=\ . NK,  r=T—\,T—2,  T—3, 


The  constraints  on  intensity  functions  for  the  backward  correspondence  are 
then  as  follows: 


B(l).  z,U)  >  0,  r«r-i,r-2 . /-I . N. 

B(2).  2,(1)  <  2,(1),  /-r-l.r-2 . /- 1,  •  •  •  .  A • 

B(3).  a,k(t)z,(t)  <  XOlk(t),  t-T-\,T-2 . /-I 


fr-1 . NK. 

B(4).  c,(M-l)r,(r)  <  ir,(T)  -  U,(t)  -  l’  c,.(r+l)- (r), 

r—r  — \ 
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t=T-2,T-3 . i=X—p+ 1, . V,  and 

c,(T):fT- 1)  ^  U,(T)  -  U'(T-l).  i=X-p-r  1 . V. 

B(5).  c,(t+ l)-(r)  <  aJtM?+\  O')  -  l'  c,(7+l)r,(r)( 

7-J  +  l 

/  =  1, . V,  all j  such  that  Dii,j)  =  1,  t=T-2,T-l, 


B(6).  a,/C,(/+l)-  (/)  ^  cap(i,l)  +  ;V,fi(/)  -  x'  o(/c((r+l)- (r), 

r-/  +  l 

/— 2 . ,V.  all  I  such  that  D(i,l)= 1,  r=T-2,T-3, 

where 


A/*(/)  =  A/w 


C-l 


«  integer  |  n  ^  X  c,(-+l)r((r) 


/=  T-X.T—2,  T—3 . 


(2.29) 


is  the  number  of  output  units  started  by  A,  after  time  t,  and 


A r,B(i)  =  A/a.v 


r-i 


//  integer |  /?<  X  c, ("+!)-,•(“) 


r-r-l.7--2.r-3 . 


(2.30) 


is  the  number  of  output  units  completed  by  A,  after  lime  t. 

The  constraints  B  ( 1 )  and  B(2)  limit  activity  intensities  to  non-negative 
values  less  than  or  equal  to  the  intensity  bounds.  The  constraints  B (3)  express 
the  exogenous  input  limitations  due  to  the  preallocations  of  resources.  The 
constraints  B (4)  require  in  cumulative  terms  that  the  output  histories  for  the 
final  products  of  the  network  do  not  exceed  the  histories 
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i=s-p+ 1 . v. 

sought  for  these  products.  The  intensity  functions  permitted  by  B<4)  need  not 
satisfy  equality  for  all  periods,  since  storage  of  final  product  is  permitted.  By 
extending  (backward)  the  number  of  time  units  used,  cumulative  requirements 
can  be  met  by  storing  final  products  produced  earlier.  Similarly,  constraints 
B(5)  require  in  cumulative  terms  that  output  histories  of  activities  producing 
intermediate  produ  do  not  exceed  requirements  demanded  by  follow-on  ac¬ 
tivities.  Again  equality  signs  need  not  apply  for  all  time  periods,  due  to  the 
possibility  of  product  storage.  Constraints  B(6)  insure  that  intermediate  pro¬ 
duct  inventories  are  not  driven  beyond  capacity. 

A  production  plan  satisfying  constraints  B(l)  through  B f 6)  is  said  to  be 
backward  feasible.  Clearly,  the  input  histories  are  not  unique.  We  shall  again 
consider  a  greedy  policy,  ie,  when  an  activity  is  called  upon  to  operate,  it  runs  at 
maximal  intensity.  Thus  the  latest  starting  times  for  activities  which  allow  the 
realization  of  the  given  output  schedule  will  be  calculated.  In  this  way  the 
shortest  time  spans  to  produce  the  given  output  amounts  are  obtained.  It  will 
be  convenient  to  keep  track  of  cumulative  activity  production  backwards 
through  time  using  the  variables 

cumfi),  /=1 . W,  r=r-2,r-3,r-4 . 

where 

t-\ 

cumfi)  =  £  c,(r+l)-  (r).  (2.31) 


The  greedy  algorithm  for  the  backward  correspondence  is  then  as  follows: 


Step  BO.  Initialize  M,B(i)  =  \B(,)  =  ami  An  =  0.  /==!, 


and  set  r=T- 1. 


Step  Bl.  For  i=\-p+ 1 . A',  set  :,(t)  = 


Min 

all  j 

such  that 
D(:.j)=  1 


For  /  =  !,  .  .  .  ,  N-p,  set  z,(t)  = 


( t ) ,  Min 

all  j 

such  that 
DU.j)-  i 


Step  B2.  Set  z^it)  =  zx(t). 


Step  B3.  For  j=2 . N,  set  Zj(r)  = 

cap(  i.j ) + N*(  i )  +a  ( /)  —  aJj  cum,  ( j ) 
aJlci(t+ 1) 

where 


|l  if  z,(r) c,(r+ 1)  +  cw«,(/)  >  j\BU) 
a(i)  =  lo  if //or  -  (2.32) 

an  indicator  variable,  shows  whether  activity  A,  completed  an  additional  output 
unit  after  time  /-I  over  that  completed  after  time  t. 


Step  B4.  Set  A B~ |  (/)  =  NB(  i)  +  a(i),  all  i, 
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where  a(i)  is  defined  by  (2.32).  Set 

(/)  =  Alfa)  +  /3(i),  alii. 

where 

1  if  :,(r)cl(r)+cumf(i)  >  1 

/3(/)  =  |o  if  not  ‘  (2-33) 

an  indicator  variable,  shows  whether  activity  A,  started  an  additional  output 
unit  after  time  r-1  over  that  started  after  time  t.  Set 

cuni'^ii)  =  cuml(i)+:l{t)cl(t+ 1),  /'=1,  .  .  .  ,  N. 

Increment  /—*/—!.  Stop  when 

cnmr(i ) 

are  all  greater  than  or  equal  to  the  desired  number  of  output  units;  otherwise, 
go  to  Step  1. 

Note  the  strong  similarity  of  the  backwards  greedy  policy  calculation  to  the 
forward  greedy  calculation,  the  former  essentially  a  "backwards  through  time” 
version  of  the  latter.  The  actual  cumulative  output  completion  schedules  for 
the  activities  are  given  by 

y,(t)  ,  ,=  T.T-\.T-2 . /-I . N. 

where 


and  p.  is  the  earliest  time  period  at  which  the  backward  greedy  algorithm  makes 
calculations.  The  system  use  of  resources  is  given  by 


/-I 

for  each  resource  A=1 . XK ,  during  each  time  period  [r,r+l). 

From  inspection  of  the  output  completion  histories 


v,U) 


t  =  T,T-\,T-2 . /  =  1 . V, 


the  backward  greedy  policy  yields  the  late  start  and  late  finish  times  for  each 
output  unit  produced  by  each  activity.  This  is  entirely  analogous  to  the  early 
start/ finish  data  generated  in  the  forward  correspondence  case. 

To  extend  the  backward  correspondence  results  further,  suppose  two 
activities  producing  intermediate  product  produce  the  same  output,  say  Aa  and 
Aa+l,  for  1  < a < N-p.  Then  the  constraints  B (5)  are  modified  for  i=a  and 
i—a+ 1  to 

B<5)\  (,)  «  a).„+1«Sl  0)  -  1'  l'  c,(r+l).-,(r). 

i—a  i**a  r=/  +  l 

all  j  such  (hat  D(ct,j)=\,  t*=T-2,T~3 .  •  •  •  . 


Note  that  coefficients  aJ  a  and  Oj  Q+1  must  be  identical.  In  order  to 
accomodate  such  substitutable  activities,  some  kind  of  policy  rule  is  needed. 
Ordinarily,  duplicating  processes  can  be  ordered  by  their  efficiency,  say  Aa+l  is 
more  efficient  or  more  preferred  for  whatever  reason.  Then  one  may  replace 
B(5V  by 


all  j  such  (hat  D(a,j)= 1,  i=T—2,T—3,  ■  ■  •  . 

With  the  addition  of  such  a  policy  rule,  a  backward  greedy  solution  is  then 
made  determinable  for  a  system  with  alternate  processes. 

The  full  spectrum  of  exogenous  input  histories  which  can  achieve  the  out¬ 
put  schedules 


l—M—p+l, 


is  only  realized  when,  in  addition  to  the  alternatives  for  preallocation  of  final 
output  and  intermediate  product  among  alternative  activities,  and  ot  exogenous 
resources  among  the  activities,  one  considers  also  the  operation  of  activities 
with  intensities  less  than  the  maximal  values  calculated  by  the  greedy  policy.  In 
particular,  reduced  intensity  assignments  allowing  more  favorable  activity  and 
resource  loading  without  violating  project  due  dates  are  investigated  in  the  fol¬ 
lowing  sections. 

In  the  special  case  that  product  inventories  are  not  capacity-constrained, 
and  intensity  bounds,  technical  coefficients,  and  resource  preallocations  are  con¬ 
stant  in  time,  the  backward  greedy  policy  is  equivalent  to  a  CPM  backward  pass 
through  the  extended  network,  with  time  assignments 
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to  each  node  (/,/??)  of  the  extended  network.  Late  finish  times  for  nodes  (/,/?:), 

i=.\-p+\ . .V,  and  the  various  values  of  m ,  are  specified  by  the  due  date 

histories 


U,{t)\,  i=N-p+\ . N. 


That  is,  a  cumulative  output  requirement  of  m  units  by  time  T  from  activity  .V 
means  node  (N,m)  has  a  late  finish  time  T.  For  each  such  finishing  node,  there 
are  critical  paths  through  the  extended  network  to  the  finishing  node  which  can 
be  determined  by  usual  CPM  procedures. 

The  above  result  also  applies  in  the  case  technical  coefficients  and  intensity 
bounds  vary  by  output  unit  (see  section  2.3),  whereupon  the  time  assignment 
to  node  (/',/«)  is  given  by 
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(2.36) 


2.5.3.  Inventories  and  The  Extended  Network 

We  have  seen  that  when  intensity  limitations  and  technical  coefficients  are 
constant  in  time,  and  there  are  no  intermediate  product  inventory  capacity  con¬ 
straints,  the  forward  and  backward  correspondences  may  be  calculated  using 
CPM  techniques  on  the  extended  network. 

It  is  also  possible  to  use  such  path  methods  wfien  there  exist  capacities  for 
such  inventories  by  making  some  modifications  to  the  extended  network.  For 
example,  suppose  activity  A ,  supplies  activity  Ay  and  the  inventory  capacity  is  n 
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units.  This  means  activity  A ,  cannot  have  completed  any  greater  number  of 
output  units  than  n  more  than  activity  A]  has  used  as  input,  ie,  node  (j,k)  must 
start  no  later  than  node  (i.k-a^+ii)  starts,  A  =  1.2 .  This  implies  all  predeces¬ 

sor  nodes  in  the  extended  network  of  O', A.)  also  precede  node  ii.k-a  +n).  If 
the  extended  network  is  modified  to  show  these  extra  precedences,  the  network 
logic  will  then  prevent  inventories  from  exceeding  capacities. 

In  summary,  for  the  dynamic  production  system  model  with  discrete  pro¬ 
duct  transfers,  the  forward  and  backward  greedy  policies  constitute  a  generaliza¬ 
tion  of  critical  path  procedures,  in  which  activities  have  variable  intensities,  out¬ 
put  streams  are  allowed,  and  product  inventories  with  capacities  are  treated. 
However,  it  shares  the  severe  limitation  with  CPM  that  a  preallocation  of 
resources  is  made,  determining  maximal  activity  output  rates. 


2.6.  Smooth  Loading 

For  either  forward  or  backward  correspondences,  an  initial  greedy  calcula¬ 
tion  will  yield  a  production  plan  with  the  shortest  time  span  to  produce  a  given 
number  of  output  units,  for  the  preallocation  of  resources  among  the  activities 
that  is  specified.  However,  such  schedules  can  have  undesirable,  jagged  activity 
loading  characteristics,  and  one  would  want  to  reduce  loading  levels  for 
smoother  behavior  without  extending  the  project  span. 

Hereafter  we  shall  assume  resource  preallocations  and  intensity  bounds  are 
constant  in  time,  and  technical  coefficients  are  also  constant  in  time  or  perhaps 
vary  by  output  unit,  indicating  learning  curve  or  other  behavior.  These  restric¬ 
tions  facilitate  the  study  of  the  smoothing  problem  in  terms  of  the  extended 
network.  Before  presenting  an  algorithm  for  smoothing,  we  shall  consider  a 
simple  example. 

Returning  to  the  simplistic  shipbuilding  network  presented  in  section  2.2, 
we  consider  the  extended  network  for  building  three  ships  with  some  specific 
technical  data  and  intensity  bounds,  as  shown  in  Figure  2.5. 

The  activity  intensity  histories  are  sketched  for  the  early  start  and  late  start 
schedules  in  Figure  2.6.  Starred  "boxes"  indicate  output  units  during  which 
activity  intensities  are  critical.  Note  that  an  activity  may  be  critical  during  pro¬ 
duction  of  certain  output  units  and  non-critical  at  other  times.  In  this  case,  the 
shortest  time  span  is  13  time  units. 

We  see  that,  in  both  schedules,  "gaps"  exist  in  the  intensity  histories, 
where  loading  ceases  for  an  interval,  and  then  recommences.  It  is  apparent  that 
serial  connection  of  activities  with  different  output  rates  causes  jagged  loading. 

In  the  spirit  of  critical  path  techniques,  one  could  consider  shifting  produc¬ 
tion  activity  between  early  and  late  start  schedules  to  achieve  less  jagged  load¬ 
ing.  Indeed,  for  this  simple  example,  it  is  evident  that  a  combination  of  early 


41 


Output  Rates:  -  cb  ~  Cc  ~  CD  ~  CE  ~  1 


Intensity  Bounds:  zA 


1  _  _  1  _  _  1  _  _  1  - 

r  Zb  ~  y  =c  ~  1'  =D  ~  2  ’  r£  " 


Figure  2.5 

Production  Data  For  Simplistic  Shipbuilding  Network 


:A  ( t )  Early  Start  Schedule  (Forward  Greedy  Policy) 


Siack: 
(A,  2)  2 
(.4,3)  3 

Slack: 

(5.1)  1 

(5.2)  1 

(5.3)  1 

Slack: 

none 

Slack: 
(D,  1),  3 
(D,  2)  2 
(A3)  1 


Slack: 
(£.  1)  3 
(£.2)  2 


Figure  2.6 

Intensity  Histories  For  Production  of  Three  Ships 
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and  late  start  schedules  can  produce  level  activity  loading  (eg,  activities  A  and  B 
early  start,  activities  C,D,  and  £  late  start).  But  such  a  schedule  does  not 
exploit  the  fact  that  activities  A ,  B,  £>,  and  £  are  noncritical  at  certain  times, 
and  therefore  resource  loading  levels  of  these  activities  could  be  reduced. 


The  duration  per  output  unit  of  activity  D  could  be  increased  from  two  to 
three  time  units  (ie,  zD  X,  :D  2,  and  :D  3  could  be  reduced  from  y  to  y) 

without  violating  the  critical  path  schedule  limitations.  A  lower  level  of 
resources  could  then  be  committed  to  activity  D  without  slipping  the  schedule. 
However,  such  a  change  would  eliminate  slack  for  activity  B .  and  thus  no 
reduction  in  intensity  of  activity  B  could  then  avoid  extending  the  project  span. 


Alternatively,  the  intensity  of  activity  B  could  first  be  reduced  from  y  to  -jy, 
but  then  the  intensity  of  activity  D  working  on  the  third  output  unit  i:D  3  =  y) 


would  become  critical. 


It  is  clear  that  there  is  no  unique  procedure  for  smoothing  unless  there  is  a 
priority  or  weighting  scheme  for  the  activities,  ie,  some  clear  mathematical 
objective.  We  shall  consider  the  natural  goal  of  smoothing  to  be  to  minimize 
activity  intensity  peaks,  where  activities  are  weighted  according  to  the  following 
formulation,  in  which  output  coefficients  have  been  generalized  to  admit  learn¬ 
ing  effects. 

Let  p,  be  the  smoothing  weight  for  activity  A,,  /  =  1 . V.  Each  p  can  be 

thought  of  as  a  capacity  cost  for  activity  Ar  Let 

« 

z\.  /-I . V, 


denote  the  peak  intensities  encountered  by  the  activities.  We  consider  a  net¬ 
work  of  A'  activities  each  producing  \M  output  units.  Using  the  notation  of 
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section  2.3,  in  which  denotes  the  intensity  assignment  to  activity  A  working 
on  output  unit  /»,  we  wish  to 

v 

Minimize  £/>,Z, 

/-l 

subject  to 

O')  zm  ^2,,  /  =  1 . Ar,  /M- 1 . V.V/. 

(//)  Z,  ^  z„  /- 1, _ M 

(///)  £  _i_  ^ 

(/./>:)  o/;  /*  ~inicim 

for  all  paths  P  from  a  starting  node  to  a  finishing  node,  where  T(P)  is  the  pro¬ 
ject  span  for  path  P  specified  by  the  initial  greedy  policy. 

(/V)  :im,  Z,  >  0,  /-i _ ,jV,  /w-i . A rM. 

The  above  program  may  be  reformulated  by  substituting  time  variables 

t,m  =  .  /=1 _ ,N,  w-1 . NM,  (2.37) 

‘-tm  C/in 

for  intensities,  and  by  making  the  substitution 

B,  =  ~r,  /«1 . /V,  (2.38) 

A 

as  follows: 

v  p 

Minimize  Y 


subject  to 
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5(1)  B,  -  cmt,m  <  0.  /-I . V.  m  —  \ . AVI/. 

5(2)  ^  4".  (=1 . V.  w  =  l, . Wl. 

5(3)  I  <  r(P). 

(/,ot)  on  P  zinicm: 

for  all  paths  P  from  a  starting  node  to  a  finishing  node,  where  T(P)  is  the  pro¬ 
ject  span  for  path  P  specified  by  the  initial  greedy  policy. 

5(4)  tjm,  B,  >  0. 

We  now  have  a  convex  objective  to  be  minimized  over  a  linear  constraint 
set.  Efficient  algorithms  exist  for  solving  problems  of  this  form,  such  as  the 
Frank-Wolfe  Linearization  Algorithm.  See  references  [1.4]. 

Many  of  the  time  assignments  t,m  may  be  critical  in  the  initial  greedy  calcu¬ 
lation,  and  so  the  smoothing  program  may  be  further  simplified  by  eliminating 
such  variables.  A  substitution  of  slack  variables 

s,m  =  tM  -  /= 1 . .V,  w-1 . N\f,  (2.39) 

Mm 

is  also  efficient,  eliminating  the  need  for  constraints  5(2).  These  steps 
effectively  reduce  the  problem  to  consideration  of  a  subnetwork  of  the 
extended  network,  namely  the  network  of  nodes  which  have  slack  for  a  greedy 
policy  production  plan.  We  shall  refer  to  such  a  network  as  the  "slack  subnet¬ 
work".  Constraints  5(3)  would  then  effectively  consider  only  slack  subpaths  of 
the  extended  network,  ie,  paths  in  the  slack  subnetwork. 

In  general,  the  program  calculates  a  minimal  set  of  activity  intensity 


% 


bounds 


4  6 


Z,U,-gr.  . v 

sufficient  to  meet  the  greedy  output  schedule.  For  our  simple  shipbuilding 
example,  only  the  peak  intensities  of  activities  B  and  D  can  be  reduced  without 
extending  the  project  span.  The  amount  of  reduction  of  each  would,  of  course, 
depend  on  the  relative  magnitude  of  the  smoothing  weights  chosen. 

For  further  improvement  in  activity  loading  for  the  given  resource  preallo¬ 
cations,  one  must  consider  extensions  of  the  project  span.  Indeed,  a  sensitivity 
analysis  could  be  conducted  by  incrementing  the  project  span  length  in  con¬ 
straints  S(3)  for  more  reduction  in  loading  peaks. 
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2.7.  Resource  Leveling 

We  have  seen  in  section  2.6  how  activity  loading  may  be  smoothed  for  a 
given  resource  preallocation.  A  smooth-loaded  production  plan  offers  technical 
improvement  over  the  greedy  plan  in  terms  of  the  loading  of  machine  services 
and  other  resources  preallocated  to  each  activity.  We  now  consider  the  prob¬ 
lem  of  improving  the  economic  efficiency  of  the  production  system  in  terms  of 
the  network  resource  loading  behavior. 

It  is  important  to  distinguish  two  different  categories  of  resources,  storable 
and  non-storable.  Storable  resources  are  materials  "used  up"  in  production, 
which  can  be  stored  from  day  to  day  before  use  at  negligible  costs.  For  such 
resources,  we  shall  assume  the  behavior  of  loading  histories  is  insignificant 
from  a  cost  standpoint,  so  that  only  the  total  amounts  used  are  of  economic 
concern.  In  considering  various  time  substitutions  of  the  greedy  production 
plan  which  meet  the  same  project  schedule,  the  integral  of  each  resource  his¬ 
tory  remains  constant,  and  so  no  economic  improvement  h  gained  from  time 
substituting  demand  for  such  resources. 

Non-storable  resources  are  the  services  per  unit  time  of  machines  and 
labor.  Such  resources  cannot  be  "hired  and  fired"  with  the  ups  and  downs  of 
system  demands,  and  consequently  resources  adequate  to  meet  peak  demands 
must  typically  be  maintained  for  the  life  of  a  construction  project.  Hence 
economic  improvement  of  a  production  plan  comes  from  time  substitutions 
which  reduce  loading  peaks  for  non-storable  resources. 

Such  motivation  for  driving  down  resource  peaks  while  attempting  to  avoid 
slippage  of  the  project  schedule  has  been  well  recognized.  Many  efforts  on 
scheduling  activities,  ie.  shifting  activities  between  early  and  late  start,  to 
improve  resource  loading  are  contained  in  the  literature.  See  references  [3.6]. 
However,  our  approach  here  will  be  to  consider  the  improvement  possible  by 


decreasing  activity  intensities  (ie.  increasing  time  assignments  to  the  production 
of  output  units).  For  a  node  (i,m)  of  the  extended  network,  the  reduction  in 
demand  level  for  resource  k  resulting  from  increasing  the  time  assignment  t.m 
is  given  by 

d  &  m:k 
im  km  ^ im 

In  the  analysis  of  resource  loading,  trade-offs  betw-een  resources  often 
arise,  w'hereby  replanning  to  mitigate  peaking  behavior  of  one  resource  may 
transfer  peaking  problems  to  another  resource.  In  the  following  analysis,  peak 
costs  for  each  non-storable  resource  are  introduced  to  properly  weigh  such 
trade-offs. 

2,7.1.  Peak  Pricing  Model 

Let  NS  ^  NK  be  the  number  of  non-storable  resources,  and  renumber  (if 
necessary)  such  resources  k  =  \,  .  .  .  ,  NS.  Let  pk  be  the  price  per  unit  capacity 
of  resource  k,  k— 1,  .  .  .  ,  NS.  Each  pk  represents  the  cost  of  maintaining  a  unit 
of  capacity  of  resource  k  (eg,  the  salary  of  a  laborer)  for  the  life  of  a  project, 
where  such  capacity  is  required  to  meet  peak  demands.  We  now  consider  the 
programming  problem  to  accomplish  a  maximal  amount  of  resource  leveling 
(for  these  peak  prices)  of  an  initial  greedy  production  plan,  without  violating 
production  due  dates.  Let 

Ak,  A--1 . NS, 

represent  the  peak  demand  per  time  unit  for  exogenous  resource  k  by  the  pro¬ 
duction  system.  For  a  production  system  in  which  each  activity  produces  NM 
output  units,  the  objective  is  to 
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vs  _ 

Minimize  £  pkAk 
k-\ 

subject  to 

1)  I  TIP ). 

( i,m )  on  P 


for  all  paths  P  from  a  starting  node  to  a  finishing  node,  where  T(P)  is  the  pro¬ 
ject  span  for  path  P  specified  by  the  initial  greedy  policy. 

R( 2)  I  I  -^8(/>,r)  <  Ak, 

w—l  /  —  I 


k—\,  .  .  .  ,  NS,  r=0,  .  .  .  ,  T—  1,  where 


h{i,m,  -) 


1  if  (/',/>;)  operates 
during  [r ,  r-1- 1 ) 

0  if  not. 


(2.41) 


*(3)  hm  >  /-I . N,  m—\ _ _  NM. 

z  c 

*/ v  mi 


R( 4)  tim,  Ak  >  0,  /  =  1 - ,N,  m= 1 . vV.V/, 


Ar-1 . NK. 


Here  constraints  R{  1)  insure  no  project  slippage,  where  the  time  assign¬ 
ment  to  node  (i,m)  is  tim  (ie,  the  intensity  assignment  to  node  (/,/»)  is  — - — ). 

t  c 

1 im v  mi 

The  constraints  R( 2)  define  the  resource  peaks  {.7;.}.  A-=l . VS.  The 

objective  minimizes  the  costs  of  resource  peaking. 

Relating  node  start  times  to  a  time  grid  r=0, 1,2 . 7"-l,  as  is  required 

by  constraints  R  (2).  is  a  difficult  task.  The  best  approach  in  the  literature  uses 
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integer  variables,  making  the  problem  too  large  to  solve  for  all  but  the  smallest 
time  grids.  (See  references  [3]  and  [6].)  The  approach  taken  here  is  to  solve  a 
sequence  of  subproblems  approximating  the  exact  formulation,  in  an  effort  to 
converge  to  an  optimum. 


Peak  Pricing  Procedure 

In  each  subproblem,  constraints  approximating  the  constraints  R( 2)  are 
applied  only  to  nodes  which  operate  inside  of  the  regions  of  peak  resource 
demand  for  the  solution  at  hand.  The  smaller,  approximate  problem  is  then 
optimized.  Inspection  is  then  made  to  see  if  peak  regions  have  shifted,  and.  if 
so,  to  redefine  the  set  of  active  constraints  for  the  next  subproblem.  This  pro¬ 
cess  continues  until  no  peak  movement  is  observed. 

Step  0.  Preallocate  resources  and  compute  the  greedy  policy  production  plan. 

Step  1.  From  the  resource  histories  of  the  current  solution,  identify  peak  time 
regions  of  each  resource  to  be  priced.  The  region  specification  is  somewhat 
arbitrary,  but  should  include  peak  and  near-peak  time  points  and  exclude  time 
regions  of  low  resource  demand.  Regions  not  priced  are  likely  to  experience 
higher  levels  of  resource  utilization  in  the  next  solution.  See  Figure  2.7  for  an 
example  of  region  pricing. 

Based  on  these  region  definitions,  compute 


8{i,m,k) 


1  if  (i,m)  operates  in  a  peak 
region  of  resource  k 
0  if  not 


/  =  ! . N,  w-1 . N\ I,  k  =  \ . NS. 


(2.42) 


If  the 


Peak  Region 


Figure  2.7 

Kxamplc  of  Peak  Region  Pricing 
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bii.rn.k) 


i—\ . N,  in  —  1 . S'M,  A  =1 , 


AS. 


have  not  changed  since  last  iteration,  stop;  otherwise.  the  subproblem  to  be 
solved  is 

vs  _ 

Minimize  £  pkAk 
k- 1 

subject  to 

(/)  £  <  W), 

(/,/>;  I  o/i  P 

for  all  paths  Pfrom  a  starting  node  to  a  finishing  node,  where  T(P)  is  the  pro¬ 
ject  span  for  path  P  specified  by  the  initial  greedy  policy. 

(//>  I  £  -2b*.  5 (/,/„,*)  4  Ak. 

/-]  Umcim 

A=1 . NS. 

( /// )  /„„  ^  ,  / - 1 . iV,  in- - 1 . AM. 

(iv)  Ak  ^  0,  /  =  1 . Ar,  /h  =  1,  .  .  .  ,NM, 

A=1 . AS. 

Equivalently,  we  may  solve 

. .  vs  \  vi/  aimk  j 

Minimize  £  £  £  pk  —  7—7  -  o(i,m,k) 

A-l  (-1  /»”i 


subject  to 
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The  subproblem  has  been  converted  to  a  nonlinear  objective  subject  to 
linear  constraints,  and  so  may  be  readily  optimized  by  algorithms  such  as 
Frank-Wolfe.  See  references  [1.4], 

Step  2.  Set  =  s„„  +  Ljm,  and  go  to  Step  1. 

At  each  iteration,  resource  peaks  are  "priced"  down,  whereby  peak  heights 
are  reduced,  but  often  peaks  are  shifted  in  location.  The  pricing  scheme 


5(/, /=! . .'V.  w-  1 . V.l/,  A  =  1 . 


is  then  changed  to  reflect  new  peak  regions  for  another  iteration.  When  peaks 
no  longer  shift,  the  procedure  terminates. 

In  each  iteration,  optimizing  the  subproblem  makes  all  nodes  in  the 
extended  network  critical.  This  can  be  seen  by  noting  that  any  slack  on  a 
start-to-finish  path  could  be  allocated  to  those  activities  operating  in  a  peak 
resource  region,  thus  decreasing  the  objective.  Thus  the  procedure  looks  only 


among  solutions  in  which  all  activities  are  critical  all  the  tune. 

The  procedure  is  similar  to  a  relaxation  method  in  optimization,  as  the 
constraints  approximating  R< 2)  are  applied  only  to  nodes  for  which  such  con¬ 
straints  are  "tight"  in  the  current  solution.  Relaxation  methods  of  optimization 
are  formally  guaranteed  to  converge  to  an  optimum  only  when  convexity  of 
constraints  is  assured.  This  is  not  the  case  here,  and  it  is  not  difficult  to  con¬ 
struct  a  simple  example  which  will  cause  the  procedure  to  cycle  two  solutions. 
However,  the  pricing  policy  may  be  readily  adapted  to  allow  the  immediate 
escape  from  such  a  cycle,  by  simply  pricing  peak  regions  from  both  solutions. 

Activity  loading  may  not  necessarily  be  smooth  from  resource  leveling. 
However,  loading  for  a  given  activity  A ,  can  be  made  smoother  by  introducing 
a  high  price  pk  on  a  real  or  fictitious  resource  k  that  is  used  by  A ,  and  not 
shared  with  other  activities.  In  the  special  case  that  each  resource  is  unique  to 
some  activity,  the  resource  leveling  problem  reduces  to  the  smooth  loading 
problem  discussed  in  section  2.6.  Further  resource  leveling  requires  a  longer 
project  span,  and  it  may  be  useful  to  study  the  trade-off  of  resource  peaks  with 
project  time  length.  The  mathematical  formulation  can  be  easily  modified  to 
address  such  a  problem  by  introducing  a  cost  rate  y  per  unit  time  on  project 
length.  For  simplicity  of  exposition,  let  us  consider  a  production  system  for 
which  the  extended  network  has  only  one  starling  node  and  one  finishing  node, 
and  for  which  the  time  T  required  to  complete  NM  output  units  is  the  length  to 
which  an  opportunity  cost  y  would  apply. 

The  Peak  Pricing  Procedure  can  be  applied  to  this  problem  by  modifying 
the  subproblem  of  Step  1  to 
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vs  .v  VI/  a,  i 

Minimize  Z  Z  Z  Pk  ~ —  7~  T  , — f  8  ( +  y  T 

k~  1  i*l  m=l 

subject  to 

.viir  £  s„„<r-  i  i„„. 

(/,/«)  o«  /*  (/,/»)  o//  /> 

for  all  paths  Pfrom  start  to  finish. 

Mil)'  s,m  ^0,  /  =  1 . N,  m= 1 . ;V.V/,  P  >  0. 

Note  that  only  one  additional  N'ariable  occurs  in  the  modified  formulation,  and 
the  remainder  of  the  procedure  could  be  applied  as  before. 

The  comments  in  section  2.6  regarding  reduction  of  the  smooth  loading 
program  to  consideration  of  only  the  slack  subnetwork  for  the  greedy  produc¬ 
tion  plan  apply  as  well  to  the  Peak  Pricing  Model.  The  subproblem  of  Step  1  of 
the  Procedure  can  be  effectively  restricted  to  nodes  which  initially  have  slack, 
as  time  and  intensity  assignments  for  other  nodes  are  fixed.  Of  course,  the 
resource  loads  for  the  entire  network  must  be  used  to  determine  peak  regions. 

2.7.2.  Application  to  Ship  Overhaul 

The  Peak  Pricing  Procedure  was  applied  to  a  network  consisting  of  151 
activities  modeling  a  ship  overhaul.  Exogenous  inputs  to  the  activities  con¬ 
sisted  of  labor  hours  of  twelve  different  craft  shops,  such  as  pipefitting,  weld¬ 
ing,  electrical  work,  etc.  For  such  inputs,  a  reduction  in  resource  peaks  was 
equivalent  to  a  reduction  in  staffing  requirements  for  the  overhaul. 

Although  the  overhaul  was  a  single  output  construction  project,  the  con¬ 
cept  of  activity  intensity  adjustment  for  resource  leveling  could  be  applied.  As 
only  a  single  ship  was  being  overhauled,  the  activity  network  would  serve  as 
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the  extended  network  for  application  of  the  model. 

Shipyard  data  on  estimated  man-hour  effort  from  each  shop  and  estimated 
durations  for  the  activities  were  available.  From  these  data,  intensity  measure¬ 
ments  and  technical  coefficients  were  developed  as  follows.  Among  the  various 
shops,  the  shop  with  the  greatest  man-hour  effort  in  a  given  activity  A  w-as 
selected  to  be  the  "lead  shop"  for  A,,  whereby  the  intensity  :ft)  of  A  on  day  / 
would  equal  the  lead  shop  man-hours  applied  on  day  i. 

The  output  coefficient  c,  for  a  given  activity  A,  was  taken  to  be 
c,  =  \\iotal  lead  shop  man-hours  needed  to  finish  A,  j  \ 

so  that 


r, (r)c,  =  {fraction  of  A,  completed  on  day  t). 

The  assumption  was  made  that  application  of  man-hours  from  other  shops 
is  proportional  to  lead  shop  effort,  so  that  input  coefficients  were  taken  to  be 
simple  ratios  of  total  activity  effort,  ie, 

(total  man-hours  of  shop  k  needed  to  finish  A,) 

a  ,  =  - - — - - i t 

(total  lead  shop  man-hours  needed  to  finish  A,) 

which  implies 

:,(t)  alk  —  (man-hours  of  shop  k  applied  to  A ,  on  day  t). 

Since  the  arcs  of  the  activity  network  represented  solely  critical  path-style 
precedences,  the  product  transfer  coefficients  were  simply  defined  by 


1  if  DdJ)3*] 
0  if  not. 


It  was  assumed  that  activity  intensities  could  not  exceed  levels  implied  by 
the  yard's  estimated  durations,  but  that  activities  could  operate  at  lower  intensi¬ 
ties,  ie,  require  longer  durations.  The  intensity  bound  tor  each  activity  A,  was 
thus  calculated  to  be 


5,  =  Cc, )  *( estimated  duration  of  A,)  1  = 

(total  lead  shop  man-liours  needed  to  finish  A  ) 

(estimated  duration  of  A,) 

As  discussed  earlier,  trade-offs  between  the  various  resources  are  a  key 
feature  of  resource  leveling,  whereby  replanning  to  reduce  peaking  behavior  in 
one  shop  may  increase  peaking  problems  in  another  shop.  To  value  the  relative 
expense  of  load  peaks,  capacity  prices  pk  were  developed.  These  prices  were 
chosen  to  be  the  sum  of  average  wage  and  overhead  rates  per  man-hour  for 
each  shop.  In  this  way,  more  incentive  was  provided  to  level-load  more 
"expensive"  shops. 

FORTRAN  IV  computer  codes  to  calculate  greedy  policies  and  to  minimize 
the  nonlinear  program  in  each  iteration  of  the  Peak  Pricing  Procedure  were 
applied  to  these  data.  An  initial  forward  greedy  policy  run  w-as  made  to  predict 
shop  loads  for  an  early  start  schedule  based  on  the  shipyard's  planned  activity- 
durations.  (Resource  preallocations  allowing  activity  operation  at  intensity 
bounds  were  assumed.)  This  run  established  a  project  time  span  of  370  work¬ 
ing  days.  The  Peak  Pricing  Procedure  could  have  been  immediately  applied  to 
the  slack  subnetwork,  but  instead,  some  observations  about  the  nature  of  the 
overhaul  served  to  reduce  the  problem  size. 

As  is  typical  of  construction  projects,  resource  loads  were  high  in  the  mid¬ 
dle  of  the  project  duration,  but  low  near  the  beginning  and  end  of  the  project. 
Activities  near  the  start  or  finish  of  the  overhaul  consisted  largely  of  tests  and 


inspections,  which,  while  requiring  time  to  perform,  were  not  labor-intensive. 
It  was  clear  that  replanning  could  not  shift  shop  workload  peaks  to  the  very 
early  or  late  stages  of  the  project,  and  that  operating  activities  during  these 
periods  at  lower  intensities  could  not  serv  e  to  reduce  workload  peaks. 

For  this  reason,  the  slack  subnetwork  under  consideration  was  reduced  to 
that  consisting  of  activities  whose  start  times  and  slack  possibly  allowed  them  to 
operate  in  the  middle  range  of  the  project  where  resource  peaks  could  occur. 
The  reduced  network  for  application  of  the  Peak  Pricing  Procedure  then  con¬ 
sisted  of  33  activities  contained  in  48  slack  paths. 

The  computer  code  used  in  the  Peak  Pricing  Procedure  applied  the  Frank- 
Wolfe  linearization  algorithm,  which  computed  a  solution  within  a  1%  tolerance 
of  the  optimum  for  each  subproblem.  (See  reference  [1].)  Peak  resource 
regions  were  defined  to  be  all  work  days  where  resource  demand  was  within  4% 
of  peak  demand. 

After  four  iterations,  peak  regions  stabilized,  so  that  no  further  improve¬ 
ment  was  possible  for  the  Procedure.  Each  iteration  had  required  about  10-15 
seconds  effective  time  on  a  CDC  6400  computer.  Figure  2.8  compares  work 
loads  of  the  mechanical  shop  for  the  greedy  production  plan  and  the  resource- 
leveled  production  plan.  As  can  be  seen,  considerable  improvement  in  shop 
loading  has  been  made.  Overall,  a  38%  improvement  in  the  objective  function 
was  made  between  the  initial  greedy  run  and  the  last  iteration  of  the  pricing 
procedure.  In  view  of  the  man-hours  involved  in  the  project  (on  the  order  of 
100.000),  such  leveling  was  significant. 

Considerable  replanning  of  activities  to  lower  intensity  levels  was  made  by 
the  Pricing  Procedure  to  obtain  this  improvement.  Undoubtedly,  there  are 
practical  lower  limits  on  activity  intensity  for  many  of  the  activities,  for  such  an 
overhaul.  However,  such  limitations  could  be  easily  incorporated  into  the 
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Time  (Workdays) 

Figure  2.8 

Mechanical  Shop  Work  Load 
Before  and  After  Resource  Leveling 


mode!  in  the  form  of  lower  bounds  (upper  bounds)  on  activity  intensity  (dura¬ 
tion).  The  constraints  .1/(1)— ,\/(2)  for  the  Peak  Pricing  Procedure  would  oniv 
require  additional  linear  constraints  expressing  upper  limits  on  slack  allocation 
to  activities  so  bounded. 

A  shipyard  ordinarily  has  a  number  of  overhaul  projects  to  be  performed  at 
a  given  time,  and  it  is  desirable  to  perform  an  integrated  planning  effort  to  con¬ 
trol  overall  resource  loading.  To  demonstrate  the  extended  network  concept, 
the  subnetwork  of  21  activities  concerning  overhaul  work  on  the  propulsion 
system  was  programmed  for  two  ship  overhauls.  In  actual  practice,  required 
overhaul  work  varies  from  ship  to  ship,  which  could  be  reflected  in  the  techni¬ 
cal  coefficients.  For  this  example,  two  identical  propulsion  system  overhauls 
were  considered,  so  that 


-/l 


and 


and 


aikl  ~  aik 2  ~  atk< 

for  each  propulsion  system  activity  A ,  and  each  resource  k\  where  r,,  c,  and  a.k 
are  the  data  used  in  the  full,  151-activity  ship  overhaul  network. 

The  Peak  Pricing  Procedure  was  applied  to  the  extended  network  of  42 
nodes  for  this  problem.  Figure  2.9  compares  workloads  for  the  shipboard 
mechanical  shop  for  the  greedy  production  plan  and  the  resource  leveled  pro¬ 
duction  plan.  This  result  is  typical  of  the  smoothing  obtainable  when  two  ships 
under  construction  interface.  Production  planning  which  is  done  merely  by 


Figure  2.9 

Shipboard  Mechanical  Shop  Work  Load  -  Two  Ships 


overlay  of  a  plan  for  a  single  ship  is  bound  to  involve  the  kind  of  unsmoothed 
peaks  show  n. 
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3.  Dynamic  Production  Planning  With  Continuous  Flow  Transfers  of  Inter¬ 
mediate  Product 

We  have  seen  how  the  Discrete  Transfer  Model  provides  considerable 
detail  concerning  discrete  interactions  between  production  activities.  However, 
providing  such  detail  necessitates  considerable  computational  complexity.  Such 
calculational  difficulties  can  be  reduced  by  considering  a  limiting  form  of  the 
Discrete  Transfer  Model,  in  which  intermediate  product  transfers  are  linearly 
proportional  to  activity  intensity.  As  will  be  shown,  the  same  form  of  linear  ine¬ 
qualities  as  in  the  Discrete  Transfer  Model  governs  this  Continuous  Flow 
Model.  The  linear  inequalities  have  been  adapted  from  those  presented  in  refer¬ 
ence  [10],  and  have  been  extended  to  include  product  inventories  and  classes 
of  exogenous  input.  It  will  also  be  shown  that  treating  product  transfers  as  con¬ 
tinuous  fiows  allows  further  development  of  the  model,  including  allocation  of 
resources  by  linear  programming. 

3,1.  Development  of  Continuous  Flow  Model 

For  the  Continuous  Flow  Model,  we  shall  initially  assume  that  each 
activity  produces  a  single  product,  to  be  shared  as  intermediate  product  input 
by  follow'-on  activities,  and/or  transferred  as  final  product.  The  case  where  an 
activity  produces  a  distinct  product  for  each  follow-on  activity,  as  with  critical 
path-style  precedences,  will  also  be  discussed.  It  will  also  be  assumed  initially 
that  no  two  activities  produce  the  same  intermediate  product.  As  in  the 
Discrete  Transfer  Model,  inventories  of  intermediate  product  awaiting  transfer 
shall  be  subject  to  storage  capacities. 

The  activity  network  will  now  be  allowed  to  contain  cycles,  representing 
more  complicated  activity  dependencies.  Leontief-like  inter-industry  flows  can 
then  be  modeled  dynamically,  as  well  as  such  phenomena  as  recycling. 
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Resources  for  exogenous  input  are  divided  into  two  categories,  storable  and 
non-storable.  Storable  resources  would  include  materials  such  as  steel,  lumber, 
fuel,  etc,  while  non-storable  resources  would  include  services  such  as  various 
kinds  of  labor  hours  or  machine  hours.  It  should  be  noted  that  limitations  in 
the  availability  of  storable  and  non-storable  resources  constrain  a  production 
system  in  different  ways. 

Non-storable  resource  limitations  constrain  the  sum  of  activity  resource 
demands  at  each  time  period.  These  limitations  constitute  capacity  constraints. 
Storable  resource  limitations,  however,  constitute  inventory  constraints, 
whereby  the  sum  of  activity  resource  demands  at  each  time  period  is  limited  by 
the  available  inventory  of  stored  resources. 

For  a  production  system  with  a  time  horizon  T,  we  shall  use  the  following 

notation  for  the  set  of  technical  limitations,  L(T ): 
w-o.1,2 . r-i 

(a)  ijit)  ,  the  activity  intensity  bounds; 

J/-1 . v 

W-0.1,2 . T- 1 

(b)  Xk(t )  ,  the  time  histories  of  non-storable  resource  exo- 

,  )A-1 . A 'S 

genous  inputs,  NS  ^  NK\ 

|r»0,  l . T- 1 

(c)  YkU)  ,  the  time  histories  of  storable  resouce  exogenous 

A-.VS+l . NK 

inputs,  where 

£r*(-> 

r-0 

is  the  cumulative  amount  of  resource  A:  supplied  during  [0, /) ; 
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(d) 


.0 


the  initial  inventories  of  activity  product  (perhaps  al 


=  1 . v 

zero)  for  intermediate  uses;  and 


(e) 


cap  (/) 


/— 0. 1 . T-] 


,  the  capacities  for  storage  of  activity  product." 


/**! . 


For  convenience,  the  intensity  r,  (/)  of  activity  An  /=1 . V,  on  each 

time  interval  [r,r+D.  / =0. 1 , 2 . 7—  1,  shall  be  partitioned  into  effort  pro¬ 

ducing  intermediate  product,  r/(r),  and  effort  producing  final  product,  -,£(/), 
where 


z]{  t )  +  ;F(  / )  =  :,(<). 


(3.1) 


Of  course,  r/(r)  or  z,F(t)  may  be  zero  for  all  t  should  activity  A,  produce  only 
final  or  intermediate  product,  respectively.  With  this  notation,  a  production 
plan 


:/(t),  :,F(t) 


/-o. 


T- 1 


/-l. 


is  said  to  be  feasible  for  L(T)  if  the  plan  satisfies  the  following  set  of  linear  ine¬ 
qualities,  A  (I  ( 7)): 

A(L(T))1.  iaik(t):,(i)  <  Xk(t),  A  =  1 . VS.  r=0. 1 . 7-1. 

/-l 

A(L(T))2.  t  I>,.(r)-(r)  ^  £  >v(r).  A-.VS-fl . \K. 

7=0  /-l  7=0 

r- 0 . 7-1. 


•  If  a  product  inventory  is  not  capacity-constrained,  the  capacity  may  he  set  euuai  to  some  \er>  large  numher 
rcr  programming  purposes 


6  6 

A(L(T) ) 21  X,  a. (“)*, ( t)  ^  inv^  +  y.  c, ( 7+1 ) r f ( ~ ) . 

r**l  i  “  1  t=0 

7  =  1 . A’,  /  =  1 . T—  1,  and 

£^(0)- (0)  <  i'hv®.  7  =  1 . A'. 

A(L(T))4.  ittVj0  +  £f,(7+l)-/(')  ~  51  1E1  ^ ( “ ) r, ( r )  ^  cap^r), 

T-0  "  r=l  ,-i 

7  =  1 . /V,  r= 1 . T—  1 ,  and 

//7V®  -  21  dl/(0):l(0)  ^  ca/?y(0),  7  =  1,  ..  .  ,N. 

/- 1 

A(L(T))5.  -/(/)  +  .-/(/)  -  r,(/),  7=1 . .V,  /=0 . 7"-l, 

Zj(t)  <  ZjU).  7  =  1 . N,  t—0 . T- 1, 

--/(/)  ^  0,  7  =  1 . /V,  /-o . r-1. 

Constraints  A(L(T))1  and  A(L(T))2  express  resource  limitations.  Con¬ 
straints  A(L(T))3  insure  adequate  intermediate  product  input  exists  to  support 
production  activity,  while  constraints  A ( L ( T ) ) 4  insure  that  inventories  of  inter¬ 
mediate  products  do  not  exceed  capacities.  Finally,  constraints  A(L(T))5  limit 
intensities  to  non-negative  values  less  than  intensity  bounds. 

In  the  case  an  activity  Aa  produces  a  distinct  product  for  each  follow-on 
activity,  the  constraint  set  A(L(T))  is  modified  as  follows.  Let 
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denote  the  inventory  of  product  produced  by  A(i  for  follow-on  activity 
.1, ,  i^a,  at  the  time  r.  Similarly,  let 

capaJU) 

denote  the  capacity  of  such  inventory.  Then  for  y=a,  constraints  A ( L ( T ) ) 3 
and  A<L(T))4  are  revised  to  the  following: 

A(L(T))3  .  iaa(r)-(r)  ^  /»r®(  +  Ica(;+1):' (:), 

r-l  r-0 

i^a,  r=l,  ....  T—  1.  and 

d,a( 0)-(0)  ^  i&a. 

A(L(T))4'.  inv®  (  +  X  cft(7+l)r£  (t)  “  X  fl,tt(7)r,0)  ^  capa  l(i), 

r- 0  r-l 

/^a,  /=1,  .  .  .  ,  T—  1,  cr/zr/ 

/>n'a./  ~  tf(a(0)r,(0)  ^  CY7Pa,(0),  /^a. 

In  the  case  two  activities  produce  the  same  intermediate  product,  say  ,4a 
and  /Ja+1,  the  constraint  set  A(LiT))  can  also  be  easily  modified.  For  this 
case  we  shall  assume  the  combined  storage  capacity 

capjt)  4-  capa^(t) 

is  available  for  storage  of  such  product.  Then  for  j=a  and  y'—a+l,  constraints 
A<L(T))3  and  A(L(T))4  are  revised  to  the  following: 
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A  ( L  ( T ) )  3  " .  ]F  X  a  Jr):  .<*)  ^  X  imf  +  £  X((:+l):'lrl, 


7=1  '  =  ] 


7-0  /> 


/  =  1 . r— 1.  and 


X °/a(o)r  (0)  ^  x 

/-I 


/-a 


/-I  o-*-l 


!  \ 


A(L(T))4”.  X^-0+X  X  c,-(r+l)r/<7)  -  X  Xa„<7>-(7)  < 


y-a 


7-0  y-u 


7=1  ,-1 


«  +  l 


X  cap.(t),  /  =  ! . T- 1.  aw/ 


X//?vy°-  X°,a^0)-,(0)  <  Zcap/O). 

y=«  /-I  y-a 


Note  that  ii  is  required  that  at  a(t)  and  a,  a+1  (r)  be  identical,  because  the 
same  product  is  produced  by  A(l  and  /4aJ.,. 

In  any  case,  A(L(T))  constitutes  a  set  of  linear  inequalities.  For  simpli¬ 
city  of  exposition,  we  shall  assume  in  what  follows  that  each  activity  produces  a 
single  product  and  that  no  two  activities  produce  the  same  intermediate  pro¬ 
duct. 

It  should  be  noted  that  in  the  case  the  activity  network  includes  cycles, 
there  must  be  non-zero  initial  inventories  of  intermediate  products  for  input  to 
the  activities  comprising  such  cycles.  Otherwise,  constraints  A(L(T))3  would 
force  zero  intensities  for  these  activities  for  all  time.  The  relationship  between 
the  magnitude  of  the  initial  inventories 


i  in 


,o 


/-!. 


and  the  final  output  realizable  from  the  network  is  discussed  in  section  3.4. 
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Although  not  required,  resources  could  be  preallocated  through  time 
among,  the  activities  as  in  the  Discrete  Transfer  Model,  and  simulation  of  the 
production  process  could  be  performed  using  greed)  policies  to  determine  the 
resulting  maximum  throughput,  or  minimum  project  span  to  reach  a  given  out¬ 
put  accumulation.  In  reference  [10],  greed)  policy  algorithms  for  the  Continu¬ 
ous  Flow  Model  were  developed  as  part  of  an  effort  to  demonstrate  the  com¬ 
plexity  of  dynamic  correspondences  and  the  dependency  upon  resource  assign¬ 
ments.  (These  algorithms  require  a  cycle-free  network.)  However,  our  approach 
here  will  be  to  consider  the  allocation  of  resources  by  linear  programming. 
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3.2.  Output  Maximization 

Assume  the  production  system  can  produce  P  ^  1  final  products.  A 
specific  product  mix  may  be  desired,  and  we  shall  consider  the  problem  of  max¬ 
imizing  the  scale  of  this  product  mix  accumulated  by  a  time  horizon  T. 

For  this  purpose,  let  rV-ri  t>e  a  variable  indexing  the  amount  of  the  pro¬ 
duct  mix  accumulated  by  time  T.  The  amounts  of  the  various  products  will  be 
related  by  coefficients 

^.V-r  1  ,/7  •  P~  ^ . P' 

where  ov+i,/>;.v+ 1  's  the  amount  of  product  p  accumulated.  These  coefficients 
specify  the  product  mix  to  be  produced.  For  convenience,  we  define  constants 

1  if  A,  produces 

Mi.p)  =  product  p  ,  /=1 . N,  p  —  1,  .  .  .  ,  P,  (3.2) 

0  if  not 

indicating  which  activities  produce  each  final  product. 

The  maximum  accumulated  output  of  a  specified  product  mix  at  time  T 
from  production  activity  during  [0,71  is  given  by  the  optimum  of  the  following 
linear  program. 

Maximize  rv+1 

subject  to 

(/)  flv+ip-v+i^  L  Z  A(/,/>)c,(/+l)r,f(/),  p=\ . P. 

/-0  (-1 

r-o . r-i 

r)  eA(UD). 

,  /-] . v 


m 


(ii)  :'(!),  :F( 
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and  rv+,  ^  0. 

In  general,  the  program  involves  (AKD  +  1  variables  (plus  slack  vari¬ 
ables),  with  up  to 

iNKHT)  +  30V)  ( 7)  +  P 

constraints.  Clearly,  the  time  horizon  (ie.  the  number  of  time  grid  periods)  is 
the  most  sensitive  factor  in  terms  of  the  problem  size  which  can  be  handled. 
The  constraint  set  A(HT )),  consisting  of  a  mixture  of  inventory  balance  equa¬ 
tions,  and  capacity  and  resource  allocation  constraints,  has  a  structure  which 
can  be  exploited  as  follows. 

Note  that  the  slack  variables  for  constraints  A(L(T))3  represent  the  inven¬ 
tories  of  intermediate  product  stored  during  the  intervals  [/,/+D. 
r=  0 . 7—1.  Calling  these  variables 

invj(t)  ,  r=0 . 7-1. 

we  may  rewrite  constraints  A(L(T))3  and  A ( L ( T ) ) 4  as  follows: 

A(L(T))3.  a,y(r)r, (/)  —  Cj(t)zj(r  —  1)  —  invj(t—l)  +  i/ivjir)  —  0, 

<-i 

j— 1, . V,  r=l,  ....  7—1,  and 

]Ta()(0)-  (0)  +  invj(Q)  =  invj*,  j—  1, . V. 

/-l 

A(L(T))4.  inxjit )  ^  capM),  y'=l,  .  .  .  ,  .V,  /=0 . 7—1. 

With  this  rewriting  of  the  constraints,  it  is  evident  that  the  constraints 
A<L(T))l.  A ( L ( T ) ) 4  and  A < L ( T ) ) 5  apply  only  time  unit  by  time  unit,  while 
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constraints  A(L(T))2  and  A(L(T))3  link  variables  through  time.  Hence  the 
constraint  matrix  exhibits  partial  block  diagonal  structure,  as  shown  in  Figure 


3.1.  As  can  be  seen, 


2(.V)(  T)  +  (A'SMD 


constraints  comprise  the  block  diagonal  structure,  with 


(AO  ( 70  +  (A K  -  NS)  (T)  +  P 


linking  constraints.  Good  potential  is  thus  offered  for  application  of  an  efficient 
large-scale  programming  procedure,  such  as  decomposition.  (See  reference 
[4].)  Using  decomposition,  at  each  iteration,  up  to  T  independent  subpro¬ 
grams,  each  with  (2Ar  4-  NS)  constraints,  would  be  solved.  In  the  case  intensity 
bounds,  resource  and  inventory  capacities,  and  technical  coefficients  are  con¬ 
stant  in  time,  these  subprograms  would  have  identical  constraint  matrices. 

As  opposed  to  maximizing  the  throughput  of  a  given  output  mix.  one 
could  consider  maximizing  the  value  of  output  produced.  We  suppose  each 
product  p  has  a  unit  price  v  which  indicates  the  value  added  in  manufacture,  or 
sales  price,  etc.  Then  the  maximum  value  of  output  accumulated  from  produc¬ 
tion  activity  during  [0,71  is  given  by  the  optimum  of  the  following  linear  pro¬ 
gram: 


Maximize  £  X  Z  vpl(i,p)c,(t+\);fr(t) 

P- 1  i-0  /-I 


subject  to 


t-o . T- 1 

z'U),  :,F(t)  e  A  (L(T)). 

J/-1 . ,v 


The  remarks  about  problem  size  and  structure  concerning  the  previous  program 


Output  Mix 
P  Constraints 


Inventory  Balance  For 
Storable  Resources  and 
Intermediate  Products 
(NK  -  NS  +  N)(T) 
Constraints 


Capacities  For 
Non-storable  Resources 
and  Inventories,  and 
Intensity  Bounds 
(NS  +  2N)(T)  Constraints 


Figure  3.1 

Structure  of  Constraint  Matrix  For  Output  Maximization 


-  n 


apply  here  as  well,  as  only  the  product  mix  variable  and  its  P  constraints  have 
been  deleted. 


Modification  for  Intensity  Lower  Bounds 

In  many  production  systems,  indivisibilities  and  other  conditions  make  the 
operation  of  certain  activities  infeasible  at  very  low  intensities,  ie,  activities 
must  either  be  "shut  off"  or  operating  above  some  minimum  intensity  level. 
While  such  constraints  are  non-linear,  the  following  formulation  converts  the 
output  maximization  problem  to  a  mixed  integer,  linear  program  for  which 
branch  and  bound  solution  procedures  exist. 

Suppose  activity  has  intensity  lower  bounds 

2,(0,  i=0, 1 . 7-1, 

perhaps  constant.  Define  a  new  binary  variable 


5 j(t)  = 


1  if  Aj  is  operating 
during  [r,/  +  l) 

0  if  not. 


(3.3) 


The  followinc  constraints  are  then  added  to  the  formulation: 


zj(t)  -  zj(t)tju)  <  o,  t-o,  i . r-i. 


Zj(t)  +  : j(t)  -  :j(t) 


1  -  6 .(/) 

J 


<  0,  r- 0 _ _  r-1. 


(3.4) 

(3.5) 


5 j{t)  <  1,  integer. 


(3.6) 


This  constraint  formulation  implies 


Sj(r)  =  0  iff  Zj(t)  =  0,  and 


For  each  activity  with  lower  bound  constraints  on  intensity,  2<  D  added  con¬ 
straints  and  T  binary  variables  are  introduced  into  the  formulation.  As  before, 
the  density  of  the  time  grid  limits  the  size  of  the  problem  which  can  be  solved 
when  such  constraints  are  present. 
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3.3.  Cost  Minimization 

In  section  3.2.  the  forward-looking  problem  of  maximizing  the  throughput 
of  production  on  a  fixed  time  interval  was  examined.  In  this  section,  the 
backward-looking  problem  of  determining  a  minimum  cost  production  plan  to 
meet  a  given  demand  schedule  for  output  will  be  investigated. 

It  will  be  convenient  to  express  output  requirements  for  each  product  type 
in  cumulative  terms.  Let 

Upit),  t= 1,2 . T,  p= 1 . P, 

denote  the  required  cumulative  production  of  product  p  by  time  t.  A  demand 
schedule  for  the  production  system  implies  constraints 

Z  X  A(/,p)c,(r-fl)-'r(7)  ^  Up(t),  p= 1.  .  .  .  ,P, 

7**0  f=*l 

r= 1 . T.  (3.7) 

Here  we  allow  final  products  to  be  produced  in  advance  of  required  schedule 
and  stored  as  inventory,  although  such  inventories  will  be  subject  to  holding 
costs,  as  discussed  below. 

Cost  Formulation 

The  objective  of  the  formulation  will  be  to  minimize  costs  for  storable  and 
m  ,  -storable  resources,  and  for  intermediate  and  final  inventories,  subject  to 
f;rr,  md  schedule  Each  of  the  cost  categories  is  discussed  below. 

•  -  ha\e  capacity  costs  corresponding  to  the  peak 

\>  me'.e  resources  cannot  be  accumulated. 
t<  icc  mmodjic  peak  loads 
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Storable  resources,  however,  have  prices.  Unless  purchasing  or  delivery 
constraints  exist,  only  the  cumulative  demand  for  such  resources  determine 
costs  for  such  resources.  If  input  coefficients  are  constant  or  vary  only  with 
cumulative  activity  output  (eg,  learning  curve  effects),  costs  for  such  resources 
are  constant  over  all  production  plans  producing  the  same  total  amounts  of  final 
products. 

Inventories  of  intermediate  product  may  -  also  have  capacity  costs 
corresponding  to  the  peak  requirements  for  storage  of  such  products.  Inven¬ 
tories  of  final  output,  however,  should  bear  an  opportunity  cost  to  properly 
discourage  holding  of  excess  inventories  of  final  products.  In  this  way,  the 
trade-off-  between  efficient  production  loading,  and  the  maintenance  of  final 
product  inventories  in  line  with  product  demand,  can  effectively  be  made. 

Assuming  linear  capacity  costs,  the  cost  minimization  problem  can  be  for¬ 
mulated  as  a  linear  program  as  follows. 

Let 

C'  -  (cf . Cfc). 

be  the  vector  of  costs  per  unit  capacity  for  non-storable  resources;  let 

c-‘-  (<:&„, . c*), 

be  the  price  vector  for  storable  resources;  let 

e:-  (c-i . c.v). 

be  the  vector  of  costs  per  unit  capacity  for  storage  of  intermediate  products; 
and  let 


Cf . CFP 
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be  the  holding  costs  for  final  products  per  unit  time. 
To  serve  as  variables  in  the  minimization,  let 


,V=  A' 


A  yc  , 


denote  the  peak  amounts  of  non-storable  resources  required  per  unit  time;  let 


^  ~  ^  .vs + 1 .  Vvj 


denote  the  total  amounts  of  storable  resources  required  by  the  project;  let 


cap  =  leap  ] . caps\, 


denote  the  required  intermediate  product  inventory  capacities;  and  let 


m\F{t)  =  (/) . invp(i)  j, 


denote  the  inventories  of  final  products  in  excess  of  the  demand  schedule  at 

time  r,  for  /  =  1 , 2 . T. 

For  given  intensity  bounds 


i=0 . T- 1 


/-I . V 


and  initial  intermediate  product  inventories 


»-i . v 


the  minimum  cost  production  plan  meeting  the  output  schedule 


/= i . r 


l/i-i . p 
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is  given  by  the  optimum  of  the  following  linear  program. 
Minimize  CM'  -f  CM  +  C:-cap  +  £  Cf-mvF(r) 

/-l 

subject  to 

Cl.  X  X  A(/,/?)c,  (r+Dc/M)  -  mvF(t)  =  L\(r), 

r-0  ,-t 

P= 1 . P,  /  =  1 . 7. 

C2.  i;^(/)-(r)  ^  A’,,  A  =  1 . A'S,  /=0. 1 . T- 1. 

/-l 

C3.  (r)  =  7A.,  A=AS+1 . VA'. 

j=*0  /-I 

C4.  £au(r)r,(r)  -  Cj(.t)ZjU— 1)  -  invjii— 1)  4-  invj(t)  =  0, 
/-l 

7=1 . ;V,  r=l,  .  .  .  ,  7—1,  and 

Xa,y(0)r,(0)  +  inyji 0)  =  /m-y\  7=1 . A'. 

/-l 

C5.  invj(t)  ^  capj(t),  7=1, . V,  r=0 . 7-1. 

C6.  z/( / )  +  =/<r)  =  :,(/),  7=1 . V.  r=0 . 7-1, 

zjit)  ^  Zj(t),  7=1 . .V,  /=0 . 7-1, 

C7.  .T  -  (.r, . ,YVS]  >  o. 

7*7 v$+i .  7  va  J  ^  0. 


so 


cap  =  t . cap\ 

^  0, 

inv'it)  =  | /'//»•{  (/),..., 

/>n  (  (r)| 

>  0. 

t=0 . 7-1, 

mvF(t)  =  |mvf  (/),..., 

mvp(l)  j 

>  0. 

/=1 . 7 

zjU),  zf(t)  >  0.  j=\ . V,  /=0 . 7-1. 


Constraints  Cl  define  the  inventories  of  final  products,  and  constraints  C2 
and  C3  define  the  peak  levels  of  resource  demand.  Constraints  C5  define  the 
required  inventory  capacities.  Constraints  C4  and  C6  deal  with  inventory  bal¬ 
ance  and  intensity  bounds  in  the  same  manner  as  the  treatment  of  the  output 
maximization  problem. 

It  should  be  noted  that  the  left  hand  side  of  constraints  C3  can  directly 
replace  Yk,  k=NS+\ . NK,  in  the  objective  function,  and  so  these  con¬ 

straints  would  not  appear  in  any  actual  programming  effort.  Also,  the  intensi¬ 
ties 


(-0 . T- 1 

/-I . v 


would  be  replaced  in  all  constraints  by  the  sum  of  intensity  effort  for  intermedi¬ 
ate  and  final  product,  ie. 


z,'(t)  +  z,FU)  =  /=1 . N,  r=0, - 7-1. 


(3.8) 


In  general,  the  program  includes  up  to 

(/>  +  AS  +  3  A')  ( 7) 

constraints.  As  before,  the  fineness  of  the  time  grid  is  the  most  sensitive  factor 


T 
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in  terms  of  the  problem  size  which  can  in  practice  be  solved.  Figure  3.2 
displays  the  structure  of  the  constraint  matrix.  As  can  be  seen,  a  bordered 
angular  configuration  is  now  exhibited.  Although  this  is  a  more  difficult  struc¬ 
ture  than  that  for  the  output  maximization  problem,  nonetheless  it  can  be 
exploited.  (See  for  example  the  discussion  of  Ritter's  partitioning  procedure  in 
reference  [4].) 

The  above  formulation  treats  a  "long  run"  problem,  in  a  sense,  because  the 
various  resource  and  inventory  capacities  are  considered  to  be  variables.  A 
"short  run"  problem,  in  which  some  of  these  capacities  were  given  data,  would 
involve  fewer  variables,  but  still  the  same  number  of  constraints.  The  degree 
of  difficulty  is  thus  roughly  equivalent  to  that  of  the  "long  run"  problem. 
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3.4.  Capacity  and  Efficiency  of  a  Production  System 

The  output  maximizing  procedure  discussed  in  section  3.2  calculates  the 
maximum  level  of  output  for  a  specific  product  mix  (or  the  maximum  value  of 
output)  producible  in  a  given  time  span  of  production.  We  can  thus  speak  of  a 
production  system  capacity  to  produce  a  specific  product  mix  (or  maximum 
value  of  products)  in  a  given  time  interval. 

The  marginal  output  return  from  increasing  a  resource  or  inventory  capa¬ 
city  can  be  studied  using  parametric  linear  programming.  (See  reference  [2].) 
Using  a  parametric  right  hand  side  on  the  output  maximization  linear  program, 
it  can  be  observed  which  inventory  or  resource  capacities,  initial  inventories,  or 
intensity  bounds,  if  increased,  would  increase  system  capacity,  and  by  how 
much.  In  general,  the  system  capacity  is  a  piecewise  linear,  concave  function  of 
each  such  technical  limitation. 

As  mentioned  in  section  3.1,  when  the  activity  network  includes  cycles, 
some  of  the  initial  inventories  of  intermediate  product  must  be  non-zero  to 
support  production  activity  by  activities  comprising  such  cycles.  The  parametric 
programming  effort  described  above  will  indicate  what  initial  inventories  are 
essential  for  system  output. 

The  optimal  tableau  of  the  output  maximizing  program  shows  which 
resource  and  inventory  capacities  are  fully  utilized,  and  which  are  not.  A 
resource  fully  utilized  up  to  capacity  can  be  said  to  be  efficiently  used,  since  a 
reduction  in  capacity  would  diminish  output  produced  in  the  given  time  span. 
On  the  other  hand,  the  tableau  shows  the  amount  of  slack  in  resource  con¬ 
straints  not  "tight".  Committing  the  full  amount  of  such  resources  to  produc¬ 
tion  activity  results  in  no  more  final  output  than  if  capacities  were  reduced  by 
the  values  of  the  slacks. 
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For  the  given  time  span  and  product  mix,  an  efficient  set  of  resource  and 
inventory  capacities  is  obtained  from  the  optimal  tableau,  namely,  the  reduction 
of  initial  capacities  by  the  value  of  the  slacks.  This  is  technical  efficiency,  in  the 
sense  that  a  reduction  in  the  availability  of  a  resource  (without  some  offsetting 
increase  in  other  resources),  would  diminish  output  producible  in  the  given 
span. 

The  cost  minimization  program,  on  the  other  hand,  determines  resource 
levels  and  inventory  capacities  which  are  economically  efficient  for  a  given 
schedule  of  final  products  to  be  delivered.  If  the  given  time  length  and  the 
optimal  output  level  from  the  output  maximization  program  are  used  as  the 
demand  schedule  constraint  for  the  cost  minimization  program,  an  interesting 
comparison  can  be  made.  The  levels  of  resource  demand  calculated  by  the 
latter  program,  which  are  economically  efficient,  could  be  contrasted  with 
resource  demands  which  w-ere  calculated  to  be  technically  efficient  by  the  former 
program.  In  general,  these  demand  levels  will  differ,  because  technically 
efficient  levels  are  dependent  on  initial  system  capacities. 

These  concepts  of  efficiency  and  capacity  are  dependent  on  the  time  hor¬ 
izon  T  allowed  for  production,  and  the  time  variability  of  input  histories 
comprising  the  set  L(T)  of  technical  limitations.  In  order  to  provide  a  connec¬ 
tion  with  steady  state  models  of  production,  we  shall  consider  the  behavior  of  a 
single-output  production  system  in  the  case  input  histories,  technical 
coefficients  and  intensity  bounds  are  constant  through  time. 

Using  the  notation  of  section  3.2,  let 

1  if  activity  A,  produces 
Mi,  1)  =  final  product 
0  if  not 


.  /  =  ! . A’. 


(3.9) 
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The  output  maximization  program  for  a  time  horizon  F  is  then 


Maximize  X  £  A (i.Dcrfit) 

i-O  /-I 


subject  to 


zjU).  zfit) 


i-0 . T- 1 


€  A  (1(D), 


J /  — 1 . V 

where  L  ( T)  consists  of  the  constant  histories 


z'\ 

t 

1*1 

H 

J 

/-l. . . 

.  ,.v 

l  J 

k- 1. . . 

.  ,  SS 

k — SS + 1 . \K 


mv 


,0 


,  and 


/-I . v 


eap, 


1,-1 . ,V 


and  the  constraints  A(L{T))  apply  according  to  constant  technical  coefficients 


k- 1 . XK 

I 

/-I . V 


j- 1. 
/-I. 


,  and 


/— l. 


Let 


.-*/(/),  r*f(r) 


/»o . r-i 


/-I . v 


be  optimal  for  the  above  linear  program.  One  way  of  looking  at  a  steady  state 
model  of  production  is  to  regard  it  as  a  statement  of  long  run  average  relation¬ 
ships.  For  such  a  model  one  defines 


{Xk,  Yk)  e  R™  -  c D(ATfc,  Yk)  6  R_ 


by 
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l'  I  U'.nc.z'fli) 

<tHXk.  Yk)  =  lim  ,~1 - - - .  (3.10) 

This  expression  is  a  meaningful  definition  of  the  long  run  average  output  rate 
•T ( A\ ,  Yk)  if  it  can  be  shown  that  the  limit  exists. 

We  consider  the  change  in  output  rate  (total  output  divided  by  time  hor¬ 
izon)  for  a  unit  increase  in  time  horizon.  Let 

i-o . T 

J /- 1 . v 

notate  the  optimal  intensities  for  the  output  maximization  program  with  time 
horizon  T+ 1.  Then 

737  £  £  -  7  rf  £  A(U>r, .-*,'<<)  - 

1  "rl  I-o  /-I  1  1-0  (-1  . 

737  £  *0. 1  +  -=t-  £ '  £  Mi.  1  )<•,-—,'■(/)  - 

1  •  1  /-I  1 ‘rl  [/-0  /-I 

l'  I  -  4=  ri  I  AUDc,.-*^/)  ^ 

f”0  /“l  /  =  ! 

737  £  M/.  l)c,.-**,f(D  +  777-7  £'  £  A</.  l)f, + 

'  ‘  1  ,«1  1  w  '  1 '  Z=0  .-1 

-r—  l‘  I  A(U)c,z**/r(r)  -  l’  i  A(/,l)c,r*ty).  (3.11) 

/  +1  i-O  /-I  i-0  /-I 

Without  loss  of  generality,  we  may  assume 

r*/( T-l)  =  0,  /—I . yV,  (3.12) 


because  these  variables  are  constrained  only  to  be  non-negative  by  AiL(T)). 


and  hence  any  optimal  solution  will  remain  optimal  if  these  variables  are  set  to 
zero.  We  claim  that 


where 


i- o. . . 
/-i. . 


.  T 
.  \ 


(.-*/(/)  for  0  <  /  <  T- 1 
(0  r=7" 


(3.13) 


and 


z'fU) 


for  0  <  t  <  T- 1 
0  r=7 


(3.14) 


is  feasible  for  the  maximization  program  with  time  horizon  7" -hi.  Since 
A(L(T+ 1))  includes  A(L(T )),  it  remains  to  investigate  constraints  for  t  =  T. 
Resource  constraints  A(L(T))1,  A(L(T))2  and  intensity  bound  constraints 
A(L(T))5  are  seen  to  be  trivially  true  for  t=T due  to  zero  intensities.  Since 


r*/(  7* —  1 )  =0,  i—\ . ;V, 


(3.15) 


inventory  constraints  A(L(T))3  and  A(L(T))4  are  also  satisfied. 
Hence 


i-0 . T 

.-1 . v 


is  indeed  feasible  for  the  program  with  horizon  7>1.  and  by  the  optimality  of 
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we  have 


I  I  AkDcr/'U)  =  l‘  I  A(/,  1  )c,z*F(i ) 
;=0  /-l  /-0  .-1 


<  I  I  A(/.  l)c,r**f(r). 

f-0  /  =  ! 


(3.16) 


Also, 


l'  I  A(/,  \  )c,z*F(t )  ^  l'  £  A ( /'.  l)c,r**/"(r)  +  I  A(/,l)c#r/, 

t=0  (=1  i=0  i  =  l  1=1 


(3.17) 


where  the  upper  bounds  on  intensities  have  replaced  /= 1 , 


is  also  the  case  that 


f=0 . T- 1 

z**f(t),  z**f(i)  ■  €  A(L(T)), 

J/  =  l . V 


since  these  constraints  are  identical  to  those  of  .-1  (1  ( 7>1 ))  up  to  time  i=T—\. 
By  optimality  of 


/=o . r-i 


/  =  !' . v 


we  have 


I  I  A(/,  \  )c,z**[(t)  ^  l’  X  A(/,  \  )c,z*fU). 
!=0  =1  1=0  ,=1 


(3. IS) 


From  (3.17)  and  (3.18).  we  have 


!l‘  I  A((,l)c,:‘f(t)  -  l'  £  A(U)c,r**fy)  <  I  A(U)c,Fr 

|/=0  /  =  !  i=0  cl  cl 


(3.19) 
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Substituting  (3.19)  into  (3.11).  we  obtain 


-r—  £  I  ^  I '  I  A</.  Uc,z*FU)\ ^ 

1  1  f-o  '1  1  1-0  =1 


1  £  a</,  i)c, r*y(  n 


7+1 


7(7+1) 


i'  £  A(/,  1  )cl:*lfr{t)  + 

i=0  i=l 


1  I  Atun-.-v 


7+1 


(3.20) 


<=1 


Replacing  intensity  variables  by  intensity  bounds  in  the  first  two  terms  on  ihi 
right  hand  side  of  (3.20),  we  have 


1  r  \  i  T—  l  \ 

ttIZ  Mi.Dc^U)  -  \  £  £  A(,.  l)c— -Of)!  < 

1  "rl  f—G  /-I  1  i-0  .  =  1  ! 


3  i  MU)c, 


7+1 


(3.21 


i  =  i 


As  7  approaches  infinity,  the  right  hand  side  of  (3.21)  approaches  zero. 
We  conclude  that  there  is  a  long  run  average  output  rate  for  a  production  sys¬ 
tem  with  exogenous  input  flows  and  technical  data  constant  in  time.  A  steady 
state  model  of  production  may  thus  be  defined  from  the  underlying  dynamic 
structure. 

In  practice,  this  suggests  an  alternative  to  the  typical  econometric  steady 
state  models  of  production,  in  which  statistical  correlations  are  made  of  system 
input  and  output  rates.  Instead,  an  activity  analysis  could  be  performed, 
whereupon  optimization  of  the  output  maximization  program  with  a  large  time 
horizon  would  provide  an  estimate  of  the  system's  maximum  output  rate  for 
given  input  rates. 
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